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JlHE following paper forms the continuation and conclusion of one on the same 
subject presented to the Royal Society last year, and printed in the Philosophical 
Transactions for 1854. I have however put it, as far as possible, in such a form as 
to be independently intelligible. 

The fourth Section (the first of this Part) contains a recapitulation of some of the 
most important results of the former Part, in the form of seven theorems, here enun- 
ciated without demonstration. 

In the fifth Section the method of the variation of elements is treated under that 
aspect which belongs to it in connexion with the general subject. It is applied, by 
way of example, to deduce the expressions for the variations of the elliptic elements 
of a planet's orbit from the results of art. 30 (Part I.), on undisturbed elliptic mo- 
tion; this example was chosen, partly because the resulting expressions are required 
in a future section, and partly for the sake of incidentally calling attention to a fal- 
lacy which has been, perhaps, often committed, and certainly seldom noticed. The 
same method, under a slightly different and possibly new point of view, is applied, as 
a second example, to the problem of the motion of a free simple pendulum, omitting 
the effect of the earth's rotation. I believe the methods of this paper might be advan- 
tageously employed in the treatment of that general form of the problem of a free 
pendulum which has been considered by Professor Hansen in his Prize Essay. I was 
unwilling, however, to attempt what might have turned out to be merely an uncon- 
scious plagiarism, without having seen the Essay in question, of which I only suc- 
ceeded in obtaining a copy on the day of writing this preface. As I now perceive 
that the investigation would be quite independent, I hope to enter upon it at some 
future time. 

The sixth Section contains some general theorems concerning the transformation 
of systems of differential equations of the form considered in this paper, by the sub- 
stitution of new variables. The most important case consists in the transformation 
from fixed to moving axes of coordinates, in dynamical problems. Some of the 
results are, I think, interesting, and perhaps new. 

The seventh and last Section contains an application of the preceding theorems, in 
connexion with the variation of elements, to the transformation of the differential 
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equations of the planetary theory. This investigation, if interesting at all, will pro- 
bably be so to the mathematician rather than to the astronomer. I think, however, 
that if the theories of physical astronomy were more frequently treated rigorously 
and symmetrically, apart from any approximate integrations; and if, when the latter 
are introduced, more care were taken to give a clear and exact view of the nature of 
the reasoning employed, it might be possible to draw the attention and secure the 
cooperation of a class of mathematicians who now may well be excused, if, after a 
slight trial, they turn from the subject in disgust, and prefer to expatiate in those 
beautiful fields of speculation which are offered to them by other branches of modern 
geometry and analysis. 

The contents of the two last Sections are more or less closely connected with the 
subjects of various memoirs by other writers, especially Professor Hansen and the 
Rev. B. Bronwin. I cannot pretend to that degree of acquaintance with them which 
would enable me to give an exact statement of the amount of novelty to be found 
in my own researches. I believe it is enough to justify me in offering them to the 
Society ; beyond this I make no claim. 

Oxford, Feb. 15, 1855. 

Section IV. 

49. The following theorems were demonstrated in the former part of this essay, 
and are recapitulated here for convenience of reference. (As before, total differen- 
tiation with respect to the independent variable t will, in general, be denoted by 
accents, which will be used for no other purpose.) 

Theorem I. — If X be a function of n variables x 19 x 2 , ...x ni and if y 19 y 29 ...y n be n 
other variables connected with the former by the n equations 

dX_ dX_ dX_ 

then will the values of x 19 x 29 ....x n9 expressed by means of these equations in terms 

of y l9 . . . . y n9 be of the form 

_dY _dY _dY 

and if p be any other quantity explicitly contained in X, then also 

*Z dY 

(the differentiation with respect to p being in each case performed only so far as p 
appears explicitly in the function). 

The value of Y is given by the equation 

Y~-(X) + (x l )y 1 +(x 2 )y 2 +... + (x n )y n9 (53.) 

where the brackets indicate that x x ...x n are supposed to be expressed in terms of 
y x ...y n (arts. 2, 3.). 
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Theorem II.— Suppose the function X to contain explicitly, besides the ^variables 
x l ...x n9 another variable t, and also n constants a 19 a 29 ...a n ; and in addition to the 
equations (50.), let the following be assumed: 

W x = bl ' m ^'da n =bn ' ••••••••• (54.) 

where 6j, ... 6 M are w other constants ; so that, by virtue of the 2n equations (50.), (54.), 
the 2n variables x l ...x n9 y l ...y n9 may be considered as functions of the 2n constants 
&i5 ... $ n > b u ... b n9 and t. Then if from the equations (50.), (54.), and their total differ- 
ential coefficients with respect to t 9 the 2n constants be eliminated, there will result 
the following 2n simultaneous differential equations of the first order; viz. — 

where Z is a function of x l ...x n ,y X9 ...y n (which will in general also contain t expli- 
citly), and is given by the equation 

In this equation -^ represents the partial differential coefficient of X taken with 

respect to t so far as t appears explicitly in the original expression for X in terms of 
x x ...x n9 a x ...a n and t; and the brackets indicate that a 19 ...a n are afterwards to be 
expressed in terms of the variables by means of the equations (50.), (arts. 5, 6.) 

TheoremllL — Let the supposition that the 2n variables x x ...x n , y x ...y n are expressed 
in terms of the 2n constants and t, be called Hypothesis I. ; and the converse suppo- 
sition that a x ...a m b x ...b n are expressed in terms of the 2n variables and t 9 Hypothe- 
sis II.; then will the following relations subsist : 



daj dyi dbj dyi 

dyi dbj dyi daj 

daj doci dbj d%{ 



!> • • { O/ . ) 



(In each of these equations the first member refers to Hyp. I., and the second to 
Hyp. II. ; and since there is no connexion between the indices of the variables and 
those of the constants, the case of i=j has no peculiarity.) 

Theorem IV. — Let the symbol [p 9 q] be an abbreviation for the expression 

^ / dp dq dp dq 
i \dy i dxi Axi dyi 

(where p 9 q are any functions of the 2n variables, which may also contain any other 
quantities explicitly ; and the differentiations are performed only so far as # 1 ,&c.,;y 1 ,&e. 
appear explicitly in p 9 q) ; then if a l9 ... a n9 b x , ... b n be expressed {Hyp. II.) in terms of 
the 2n variables and t, the following equations subsist identically : 

[a i9 ft,] = - [b i9 aj = 1, [a i9 ft y ] = [a,, ff y ] = [b i9 ft,] =0 • . . . (58.) 
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(i being different fvomj); and obviously in all cases 

[p> q]=-[9>p]> and Ip>p]=° ( art - 9 -)- 

TheoremN. — If u 9 v be either (1) any two functions whatever of the 2n constants 
a 19 &c, b 19 &c, or (2) any two functions whatever of the 2n variables x 19 &c, y }9 &c. 
(which may in either case also contain t explicitly)*, then 

^ \ du dv du dv 



dyi dxi d%i dyi 



du dv du dv i (zq \ 

da* dbi dh da^* * * * * * * * V •/ 



(When u, v represent functions of the constants, the differential coefficients in the 
first member of this equation refer to Hyp. II.; and, when functions of the variables, 
those in the second member refer to Hyp. I.) (art. 10.). 

Theorem VI. — Let x l9 ... x m y l9 ... y n be 2n variables concerning which no supposi- 
tion is made except that they are connected by n equations of the form 

a i z=z(p i yx l9 x 2 , .... x n9 y l9 y 2 ^ ••• y<n) (a.j 

(where the functions on the right are only subject to the condition that the n equa- 
tions (a.) shall be algebraically sufficient to determine^, ...y n in terms of x l9 ...x n9 a 19 
&c, and may contain explicitly any other quantities besides x 19 &c, y l9 &c). 

Then, if by means of the equations (a.) the n variables y 19 y 29 ... y n be expressed as 

functions of x 19 x 29 &c, a 19 &c. ; in order that the — ^ — conditions 

dyi dyj 

1X00 j (X/Uu^ 

may subsist identically, it is necessary and sufficient that each of the -^ — - expres- 
sions [a i9 dj"] vanish identically. 

Theorem VII. — Let Z be any function whatever of 2n variables <x x ... x n9 y x ...y n9 and 
i. If of the system of 2n simultaneous differential equations of the first order 

, dZ t dZ /T N 

*'=*<'*= "a; < L > 

there be given n integrals involving n arbitrary constants a 19 a 29 ... a n , so that each of 
these constants may be expressed as a function pf the variables x 19 &c, y 19 &c. (with 

n (n — 1 ) 

or without t) ; then if the — ^ — conditions [a i5 aj =0 subsist identically , the remaining 

n integrals may be found 9 as follows. By means of the n given integrals let the n 
variables y x ...y n be expressed in terms of x l9 &c, a l9 &c. ; and let (Z) be what Z 
becomes when y x ...y n are thus expressed. These values oly l9 y 2 ...# w and — (Z), will 
be the partial differential coefficients with respect to x l9 x 29 ... x n and t 9 of one and the 
same function ; call this function X, then, since its partial differential coefficients are 

* It was inadvertently stated in art. 10, that u, v must not contain t explicitly. But it is evident that no 
such limitation is implied in the demonstration of the theorem. The preceding theorem is obviously a particular 
case of this ; namely, the case in which u=aj, v=bj. 
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all given (by the equations -^=y i9 ~^-=-- (Z) J, X may be found by simple inte- 
gration, and is therefore to be considered a given function of x 19 ... x n , a 19 ... a n and t. 
The remaining n integrals are then given by the n equations 

b x ... b n being n new arbitrary constants. 

[On the relation between this theorem and the theories of Sir W. R. Hamilton 
and Jacobi, see arts. 15-20.] 

50. Other results established in the former part will be referred to as occasion may 
require. To the theorems enunciated in the preceding article, the following may 
now be added. 

Returning to the equations (50.), (54.), (55.), we may observe, that if, in thejirst 
members of (55.), x { , y { be supposed expressed in terms of a l9 &c., b l9 &c. and t 9 then 

7 7 

— , -| £ may be written instead of x i9 y\ ; since on this hypothesis the total differential 

coefficients of x i9 y { are obtained by differentiating with respect to t as it appears 
explicitly. We have therefore 



IXiAJ ' 



dZ dyi _^Z 



dt dyi dt dxi 

where the first members refer to Hyp. I., and the second to Hyp. II. But since the 
equations (50.), (54.) involve a 9 b 9 exactly in the same way as they involve x 9 y, it is 
obvious that the same reasoning which leads to the equations just written, would 
lead, mutatis mutandis, to the following, which may be considered as an addition to 
the system of equations (57.) (Theorem III.): 

da* dZ dbi dZ 



dt dbi dt dcti 



(60.) 



In these equations, a i9 b ( in the first members are supposed to be expressed in terms 
of the variables {Hyp. II.), whilst in the second members x 19 &c, y 19 &c. are supposed 

dH. 
to be expressed in terms of the constants and t (Hyp. I.). As before, Z= — -^p but 

in (60.) Z is differently expressed, being what the Z of (55.) becomes when x u &c a , 
y l9 &c. are expressed according to Hyp. I. 

It is to be remembered that all consequences deduced from the form of the system 
(50.), (54.) belong to the system of equations, obtained as in Theorem VII., which 
express the solution of the differential equations (I.). Such a solution will be called, 
as before, a normal solution; and the system of equations obtained by expressing 
a l9 &c, & 1? &c. in terms of the variables and t 9 will be called a system of normal inte- 
grals. (See art. 20, and the note to art. 29.) 

51. Let a l9 &c, b 19 &c. be called, as before, elements. If then c be any function of 
the elements, when the latter are expressed in terms of the variables and t 
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{Hyp. II.), c becomes also a function of the same ; and we have 

dc y / dc ddi^dc db\ y (dc dZ^dc dZ\ fai \ 

dt \da,i dt dbi dt J \dcti dbi dbi daj 

(see the last article). But, by Theorem V., this becomes 

ihM- • • • • • • («a.) 

It is worth observing that both this equation and (60.) might have been obtained 

dc 

indirectly as follows. Since c is constant, we have c'=0 ; that is, -g+[Z, c]=0 (see 

(32.), art. 22.); this gives (62.), since [Z, c] = — [c, Z], and again, by Theorem V., is 
changed into (61.); and if, in the latter, we put successively c=a j9 c=^, we obtain 
the system (60.). 

Section V. — On the Variation of Elements. 

52. The following general problem includes, I believe, all the cases which occur in 
practice. Let P 13 ... P w , Q 13 .... Q n be any functions whatever of the 2n variables 
x 19 ... x M y 19 ... y n and t. It is required to express the 2n integrals of the system of 
2n simultaneous differential equations of the first order 

*;=?,, ^=0, (63.) 

in the same form as the integrals (supposed given) of the canonical system 

/ dZi i d/i /t \ 

^=-7-, ^= — — ? (1.) 

1 diji ^ l dXi v ■ / 

by substituting functions of t for the constant elements of the latter system. 

Suppose a normal solution (see end of art. 50.) of the system (I.) to be employed. 
The elements a^ b { represent the same functions of x 19 &c, y 19 &c. and t as before, but 
are now variable ; consequently we have 



/ dcti I ddi , dCL{ , 

a s - 



dt +i 4^ y ^j — *" 4 ^ y ^ 



with a similar expression for £-. But, by equations (57.) an d (60.), these are imme- 
diately transformed into the following : 



b'. = — — — 2.{q.^— P-^ 



(E.) 



where Z, Q,, P,, #,, y s in the second members are supposed to be expressed (Hyp. I.) 
in terms of the elements and t. Thus the system (63.) is transformed into a system 
involving the new variables a i9 b i9 instead of the original variables x i9 y t . 

53. If, instead of employing a set of normal integrals of the pattern system (I.), we 
take any complete set of integrals c 19 c 2i ... c 2n , then c l9 &c. may be considered as 
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functions of a u &c, and again, through them, of the variables. We have then 

c *—da l ai ^'"^db 1 ° 1 'T m --- ; 

and if in this equation the values of a' l9 &c. be introduced from the formula (E.) of 
the last article, the following expression results : 

c;={Z, c*} +2j(Qj{x J9 cj — T?j{y j9 cj) 
(in which the symbol {p, q] is used to denote 

y \ dp dq dp dq\ 
\db k da k dak db/cj 

so that by (59.) (Theorem V.) we have {p, q} = — \_p> g] ; but in {p, q} p and q are 
considered as functions of a 19 &c, b 19 &c, whilst in [p 9 q] they are considered as 
functions of x l9 &c., y l9 &c). Now, considering^, q as functions of c 19 &c, and through 
these, of a 19 &c, we have (by reasoning exactly similar to that employed in deducing 
equation (24.), art. 9.) 



<»*H { ^ J GU-££)) 



(the summation referring to all binary combinations of the indices a, j8). Hence we 
have, putting q~c i9 

{P,c i } = 2j{c a ,c i }£-J, (64.) 

and consequently the above expression for c- becomes 

c t i ^{Z 9 c i }+^^{c a9 *}(<%$-?$)), ....... (F.) 

an equation which is easily seen to become identical with (E.), art. 52, when c x ...c 2n 
represent a x *..a M b x ...b n . 

54. The simplest case is that in which the system of equations (63.), whose inte- 
grals are sought, are of the canonical form ; that is, where 

dyi l dxi 

W being a given function of the variables (with or without t). In this case the 
formula (E.) becomes 



1 dbi dbi 
h ,_ dZ.dW 



> 



(65.) 



dcii da,} - 

whilst (P.) is easily found to be reducible, by the help of (64.), to either of the fol- 
lowing forms : 

c;={Z,c 4 }-{W,c,}. (66.) 

c;=x({c a5 cj(g^)) . . .. (67.) 
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If we put W=Z+0, so that £1 may be called the " disturbing function," the above 
formulae become 

/y'_. ®^ "h! ff\Q \ 

%J,t m " — — — ^ f / • ■ . . a , » » • a a • s a I vl O • / 

c\-=^ a ( {c i9 cj^r- 1 ....... (by.) 

On the first of these forms see the note to art. 38. With respect to the form (69.), 
if we put for {c i9 cj its equivalent — [c„ cj, or [c a5 cj (see Theorem V. art. 49.), we 
obtain the well-known expression 






The difference between this last form and (69.) consists in this ; that in the latter the 
coefficients [c a9 c 4 ] are obtained from the expressions for c l9 c 2 , &c. in terms of the 
variables; whereas in (69.) the coefficients {c i9 c a } are similarly obtained from the 
expressions for c 1? &c. in terms of the normal elements a l9 &c, b 19 &c.* ; and when a 
normal solution of the undisturbed problem has been obtained, the latter process will 
generally be found much more convenient than the former, since the elements c 19 &c. 
will usually be much simpler functions of the normal elements than of the variables. 

55. In illustration of this, it will be worth while to deduce the expressions for the 
variations of the ordinary elliptic elements of a planet's orbit from those of the normal 
elements given in art. 30. 

Let a and e be the semiaxis major and excentricity, ; the inclination of the orbit to 
a fixed ecliptic, v the longitude of the node, w the longitude of the perihelion, w#+(e) 
the mean longitude of the planet ; longitudes being reckoned in the plane of the 
ecliptic (from a fixed origin) as far as the node 9 and then on the plane of the orbit. As 

usual, n stands for~r- Also let ^+(g) = l ndt-\-$, so that s f =(s) f +tn f . 
If, then, we call the six normal elements k 19 co 29 a 3 , j3 w |3 3 , /3 35 we have (see art. 30.) 

mp r> _ / sy— (s) 



-i=^» ft 



2a 



n 



a 2 =ms/pa(\ — e 2 ) 9 /3 2 :=sr— -i>, 

«3 = ^ N / 1 M(r^).COS/, /3 3 =y; 

from which, conversely, 

1 e —mV Ps? 



«o , , _ _ (2« 2 )^ 



; 3 






* 



\abj daj daj dbjj 
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From these expressions the values of {a, e}, [a, i} 9 &e. are found with the greatest 
simplicity, and the results are 



mp [a, (g) } = 2 na 2 , mp { (s) , e } 



na s/ \ 



,2 



(1-x/l-e 2 ), 



mp{w y e} =— * m(jj{(t) 3 /}=— 7=== tan~> 



Vl — e 



,2 



m(A{&,t} =^r|=itani 5 m^{v) 



^« 



in * V I — e 



sin 



,2' 



the rest all vanishing. Hence, observing that if R be taken in its usual signification 
e have = — R, we obtain # 



pa! 



(Jbe 



a(e) 



{S+(l-s/W)f} 



M0' : 



2wa 2 ^+ 

da 



dR . na^/i-i- e 2 



^R 



, dR , B£ , 1 //^ 






tan 



* rfRl 



e rfe Vl — <? 2V "*2 



5T/' 



n« 



f =VT-^ 



na 



p 



_1. ^5+ tan if— +— 

sin 1 dv ■ 2\d(s) dm 

dR 



>j 



JR 



sin * v^ 1 — e 2 di 

in which we may, as usual, put s for (s), provided that in forming the term ^ 3 nt be 

exempt from differentiation with respect to a. 

56. A comparison of the above process with that by which the corresponding 



* If we consider R as a function of p, q instead of t, y, where _p= tan t cos y, g=tanf sin y, we find 

dR 



dt 
dR 



= see 2 * f 



dR . dR\ 

cosy-T-+smy-3rr) 



dp 



dq 

dR 



and consequently 



I . ( dR . dR\ 
— = tan * ( cos y — sin v — ]» 

dv \ dq dp J 



, —na(sect) q \ dR , , t 
y*p = — - v ■■{ sec t -- r ~ + tan - cos v 



Vl—e* 



dq 



\d(s) dm) 



,_rca(sec*) 2 f dR 1 . (dR dR\ 

w -vr^{ sec ^- tan 2 smv v^ + Wj* 

The formula will then agree with those of the Meeanique Celeste (Supplement to vol. iii. p. 360, ed, 1844), if 
we allow for the different mode of measuring longitudes, and neglect, as Laplace does, terms of the second 

order with respect to 1 and f?5. (Laplace uses R with the opposite sign.) Those in the text agree (allowing 

for notation) with the expressions given by Professor Hansen, Astr, Nachr. No. 166, art. 3, equations (2) 
MDCCCLV. 2 T 
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expressions are obtained by Pontecoulant*, will show the convenience of using the 
coefficients {c i9 Cj} instead of [c i9 cj] (in Pontecoulant's notation (c {> c y )). 

[It will be observed that the formulae for (s)', vf 9 i at the end of the last article, 
do not agree with those of Pontecoulant (p. 330) for the variations of the corre- 
sponding quantities g, a, <p. The reason of this is as follows : — In Pontecoulant's 
notation <p expresses the same as 1 in this paper, and a the same as v. But &> (the 
longitude of the perihelion) is not the same as «r; the former being measured entirely 
in the plane of the orbit from a radius vector, fixed in that plane\, and assumed as 
the origin of longitudes. Consequently g, in Pontecoulant (which we will call g, for 
distinction), is not the same as (g) in the present paper. In fact, if we equate the 
expressions for the mean anomaly in the two notations, we have 

also it is evident that if we put /3 for the angle between the node and the origin from 
which & is measured, we have dfi= — costdv, and nr=i>+j3+0, so that 

d&=dco-\- (1 — cos $)dv. 

If then it were allowable to consider R as capable of being expressed as a function of 
& and g, instead of m and (g), and if we represented by (R) the expression for R so 
transformed, we should have 

and if, in the two first terms, we put for dW and rf(g) the values £for=rffi;+(l— cos*)A, 
rf(g)=rfg i +(l— cos/)cfr, and compare the two expressions, we find 

^R rf(R} ^R rf(R) 

^R / X /^R dR\ rf(R) 

— - I- ( 1 — COS l ) ( "77T -4- -j — I = -—, — • 

These relations, together with the equation 

fflr , = fl/+(l— cos/)^, 

are easily seen to render the expressions at the end of art. 55 identical with those of 
Pontecoulant ; in fact, it is by an equivalent transformation that the latter are 
finally obtained by that author from the correct expressions in p. 328. But it is to 
be observed that this proceeding is founded upon a false assumption; for it is im- 
possible to express R as a function of a, e, /, v, g„ a>, as is obvious from the considera- 
tion that R, in its original form, is not a function of (g)— m merely, but also of (g) ; 
whilst (g) is not expressible as a function of the new elements, as is shown by the equa* 
tioni(g)==rfs ; +(l— cos i)d»$. It would be out of place to enter further into this sub- 

* TMorie Anal, du Systeme du Monde, tome i. pp. 316-330. 

f On the meaning of this expression, see below, art. 73. 

I It would be a work of some trouble to trace accurately the process by which Laplace arrives at the for- 
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ject here, especially as some of the most important principles involved in it have 
been discussed elsewhere # . See also Appendix B.] 

57. Returning to the expression (69.), art. 54, it may be observed that the coeffi- 
cients {c i9 Cj} are to be expressed in terms of c l9 c 2? &c., and this involves no difficulty 
when each of the two sets of elements c l9 &c., a 19 &c. can be expressed in terms of the 
other explicitly, as was the case in the example just discussed. Suppose, however, 
that the normal set a 19 &c, b 19 &c. are given in terms of the set c 15 &c, but that it is 
impracticable or inconvenient to obtain the converse equations expressing the latter 
in terms of the former. In this case we may proceed as follows. 

Adopting the notation of art. 1-J-, and putting/, g for any two of the set c l9 c 2 , &c, 
we have 

suppose this equation written at length, and then, after multiplying by j h a - j , let 

each side be summed with respect to all binary combinations/,^. The result is (see 
art. 1, equation (4.)), 

2(^-{/,^}) = l • (70.) 

(the summation referring to the combinations / g). Again, if the former equation 

d{ 7) 0) 

be multiplied by Yl ¥ , where p 9 q represent any two of the normal elements, a l9 &c. 
b l9 &c, except a conjugate pair, and the sum be taken as before, we have 

2 (w {/ ^ =0 (7L) 

The two formulae (70.), (71.) give n(2n— 1) linear equations for determining the 
n(2n—l) unknown quantities {f 9 g}; the coefficients of the latter being all given 
functions of c l9 &c. But such cases will hardly occur in practice. (With respect to 
the form of the above system of linear equations, it is easy to show that the complete 
determinant of the coefficients is =1.) 

58. The integration of the formulae (65.), art. 53, would give the means of ex- 
pressing the solution of the system 

,_ dW /__^ 
i- "~ diji ? y % ~ "~" dxi 

mulse alluded to in a preceding note, as the various steps of it are to be found in different places, the notation 

is somewhat inconsistent, and the results do not profess to be rigorous. My impression is, however, that 

Laplace nowhere commits the fallacy of assuming (for example) that R is a function of r, v, z, or r, v, s (see 

vol. i. p. 295), where v is the angle described by the radius vector on the varying plane of the orbit. 

* See Jacobi's two letters to Professor Hansen in Crelle's Journal, vol. xlii. 

, . . d(u, v) , , • J .. r du dv du dv 

f 1. e. using -) ■<- as an abbreviation for -— -— 

d(oc, y) dx dy dy dx 

2 t2 
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in the form of a normal solution of any other similar system 

/ dZ , dZ 

1 dyl y i dxi 

which may be chosen as the pattern. 

In the most usual examples the function to be chosen for Z is naturally suggested 
by the circumstance, that W presents itself under the form of the sum of two func- 
tions Z+f2, of which the former, taken alone, gives an integrable system. But this 
is not necessarily the case; and it is worth while to observe that the forrnulse (65.) 
take a simple and remarkable form whatever Z may be, provided that it be a func- 
tion not containing t explicitly. For then, assuming the "integral of vis viva," Z=k f 
as one of the normal integrals of the pattern system # , the element conjugate to h is 
r (the constant added to t) ; and observing that Z, in (65.), being expressed in terms 

dZ 
of the elements, reduces itself simply to A, we shall have ^=1, whilst the differential 

coefficients of Z with respect to all the other elements vanish ; so that, if we put 
«w*«n-» ^•••^n-i for the remaining elements, the system (65.) takes the following 

f ° rm :_ u rfW , , , dW 

n — dr' 7 1+ dh 



' dbi * das 



<"• » "» • » • . • [/ 4U , I 



This, in dynamics, gives the process to be used in the following problem : " To express 
the solution of any dynamical problem in the form of the solution of any other (involving 
the same number of variables) in which the principle of vis viva subsists." 

59. As an example of the above process we may apply it to determine the motion 
of a simple free pendulum (not taking into account the earth's rotation). 

Let / be the length of the pendulum, and let the mass of the material point m 
placed at its extremity be represented by unity. Also let x 9 y, z be the rectangular 
coordinates of m, the origin being at the position of rest of m, and the axis of % 
directed vertically upwards. The equation to the sphere described by m is 

x 2 +y 2 -\-z 2 —> 2fe = 0, 

and the force function U is —gz. 

Hence if we take, as the two independent coordinates, the radius vector g of the 
projection of m on the plane of xy, and the angle 6 between g and the axis of «r, we 
shall have for the differential equations of motion, 



du dv 



i dW , dW 
da dv 



> 



(A.) 



* See art. 19 (where h { in equation (29,) is a misprint for b;). 
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where u, v, are the variables conjugate respectively to §, 6, and defined by the equations 



u 



dT 






dT. 



df " dff 
and W is T— U expressed in terms off, 6, u, v. 



Now x=g cos d, y=g sin 6, z=l—\/l*—f ; hence 



T[=l(x' 2 +y' 2 +z 



*)Ui(— 



g 



,f+g 2 <? 2 



from which the following expression for W is easily obtained : 



w=^ 2 - ? 



2\ I 



2 



*-u*+ v ^)+g(i-Vi 2 -?). ...... (W.) 



Now let us take as a model for the solution of the above system, a set of normal 
integrals (in polar coordinates) of the system 

x"+n 2 x~0, y n +n 2 y=:0, (B.) 

where n 2 =f. In this system we have U= — -$n 2 f\ and proceeding exactly as in 
art. 27, we obtain the following results : the two integrals of vis viva and of areas are 



,2 



h~v(u 2 +^+nY 



2 
c=v 



0-) 



these are to be solved for u, v\ and then V is to be obtained from the equation 
dV =udg-\-vdd. This gives 



V=c0+fa 



2h-ny -*)*; 
§ J 



and the remaining integrals are given by the equations 

dY . . dV 
Th =t+T > Tc =7 *> 

r and m being the elements conjugate respectively to h and c. Performing the dif- 
ferentiations first, and taking the integrals in the second term so as to vanish with 

the expression \2h~- n 2 f— %Y (see Appendix A.), we find easily the final equations 

n 2 f=h+*/tf-- rcV.cos 2n(t+r)} 
&f 2 =.h—s/h 2 — n 2 c 2 . cos 2(d~«?)j 9 

in which w is the angle between the axis of x and a (distant) apse, and — r is the 
time of passage through that apse. The four equations (i.), (ii.) comprise a com- 
plete normal solution of the equations (B.). The last is the polar equation to the 
elliptic orbit ; and if we call a, b the semiaxes of the ellipse, we have 



(ii.) 



c=nab, h=n: 



2 
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60. The solution of the system (A.) of the last article will now be expressed by the 
same equations (i.), (ii.), if the elements h, c, r, & be variables defined by the system 

(see art* 58.) 

n dr> T ~ 1 + dh 

( j = _dW m i = dW 
dvr f dc 

where W is to be obtained by substituting* in the expression (W.), art. 59, the 
values off, 9 n> v in terms of the elements and t, derived from equations (i.), (ii.). 
The result of this substitution is 



+^WT C0S 4n{i+r) +nV ( / - ^ /F= ?)^ 

in which the value off 2 in the last term must be understood to be substituted from 
the first of equations (ii.). if we call <p the angle between the pendulum and the 
vertical, we shall have evidently 

f?l(l-y/K^f) =»V( 1 - COS <f>), 

and the differential coefficient of this term with respect to any constant k involved 
in the value off will be ~- ~ ^- Observing this, we obtain the following expres- 
sions for the variations of the elements : 

ti=—n(Mcp—l)s/7iF^ 

Ml 

h 






+ I(sec ? -l)(l+- 7F i=co S 2»(^+r))+-A ? cos4^+r) 



c' = 

^=^--(sec <P-l)~^~^ cos 2n(t+r) -^ cos 4n(/+r). 

The third of these equations gives ab= constant ; hence, by means of the equations 
at the end of art. 59, the following expressions are easily deduced : 

-= — - j=g j — ( sec ^~" ! ) sin 2n ( t +*)+—]*-' sin 4n(t+r) 

% I TO 1 > Q ■ JO \ 

r '=^r-(~~ ! + cos M^+r))-f-(sec 0- 1 ) U +^z^cos2;?^+r) J 



Jzzmabl- ^ *— (seep— 1) — ^jg— J- 



These equations are rigorous, and in general not easier to integrate than the original 
system of which they are a transformation ; but they may be integrated approxi- 
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mately on particular hypotheses. For instance, if the pendulum never deviates much 
from the vertical, ^ is always small, and sec <p— 1 = |^ 2 nearly ; introducing this value, 
and substituting for f, we have for the non-periodic parts of the above expressions 

i a? + b 2 ( 3 nab 

16P 8 l l 

Hence the axes of the mean ellipse are constant; also if T be the time of describing 

2tt 1 

the ellipse, we shall have approximately T=~-j— -,-, whence it follows that the mo- 
tion of the apse during this period will be -|^-- * r > ov -^-2tt nearly. 

This agrees with the statement of the Astronomer Royals The above approxima- 
tions would cease to be sufficiently accurate if b were very nearly equal to a, or the 
motion very nearly circular ; but they apply to any other ease, and in particular to 
that in which b is very small, or the motion nearly rectilinear. In the case of nearly 
circular motion, it would be necessary to develope sec <p in a series of powers of 
a 2 —b 2 , and the results would be applicable whether <p were small or not. But I shall 
not pursue this subject further here. 

Section VI. — On the Transformation of Variables. 

61. The method of the variation of elements, theoretically considered, consists 
merely in a transformation of variables of a particular kind ; that kind namely, 
which leads to a new system of differential equations belonging to the same general 
class as the original system. But practically, the choice of variables is determined 
by the well-known considerations from which the method derives its name. 

It is the object of the present section to consider the general class of trans- 
formations of which the method in question is a particular, and not the only useful 
case. 

62. Definition of Normal Transformations. 

Let |j, | 2 , ... | w , fi l9 fi 29 ... vj n be new variables connected with the original variables 
x l9 &c. y u &c. by 2n equations (which may also involve t explicitly), such that each 
variable of either set may be considered as a function of the variables of the other 
set (with or without t). Let P be any function of §„ % 29 ... \. m y l9 y 29 *..y n and t ; then 
if the equations connecting the old and new variables can be put in the form 

■ VL XT (X/jl * ^ i 

I propose to call the transformation normal. 

* Proceedings of the Royal Astronomical Society, vol. xi. p. 160. 
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[It follows from Theorem I. art. 49, that the system (73.) is equivalent to each of 

the following: 

dQt dGt 

Ifa—y* djf i :=: ~ m9li > 

in which Q= — P+2j(#<y*)> and is' expressed in terms of x l9 ... o? n , g 2 , ... £»; or 

^R „ ^R 

in which R= — P+2i(5#«), and is expressed in terms of y l9 ... y n9 q 19 ... n n \ or lastly, 

in which S= — P+S^^+l^), and is expressed in terms of x 19 ... x n9 vi l9 ... yj n . 

Any one of these forms might be used ; but I shall employ the form (73.) for 
reasons of convenience.] 

63. Inasmuch as the equations (73.) of the last article are of the same general 
form as the system (54.), art. 49, all the conclusions deduced from that form will 
subsist, mutatis mutandis ; so that we may apply the Theorems (II.), (III.), (IV.), (V.), 
art. 49, by merely changing X into P, and 

x l9 ... x n9 y 19 ... y n9 a l9 ... a n ,b„ ... b n9 respectively into 
observing that instead of x i9 y t we must now write -g> -r- l# - We thus obtain the fol- 

at dt 
lowing relations : 

dh_dW- d^_ _d^ ffJA , 

dt~~ d* dt~~ &$i •••••••■ •■ • (74. ) 

where Y is a function of l l9 &c, ti l9 &c. and t 9 defined by the equation 

T — - I ~17 J } • • • » ....... \jD.J 

the brackets indicating that the expressions for y l9 ... y n m terms of the new variables 

dV 
| x , &c, yj 19 &c, are to be substituted in -^ after the differentiation with respect to t; 

which is performed so far as t appears explicitly in the original expression for P as a 
function of \ x ... l n9 y Y ... y n and if. (See Theorem IL) 
We have also the svstem 

d$i dxj dii^dyf 



dyj . dt\i dxj df\i 
dvji __ datj dyi dyj 



?« 



dyj d%i dx § dgi 



(76.) 



For In the original theorems x { is the same thing as the differential coefficient of w% taken with respect to t f 
as t appears explicitly in the expression for w\ in terms of a v &c., b v &c. and t; the analogous quantity in the 
present case is therefore the differential coefficient of £ taken with respect to t t as £ appears explicitly in the 
expression for & in terms of # p &c, y p &c. and £. But this must not now be denoted by ^, inasmuch as 
x l9 &c, y x > &c. are themselves afterwards to be considered as functions of t. 
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(see Theorem III.). And if p, q be any two functions of the variables (with or with- 
out £), then 

\dyi dxi dxidyj \dvji dfa d^drn)* ' ' " ' U/v 

where p and q in the first member are supposed to be expressed in terms of x l9 &e., 
y u &c, and, in the second, in terms of f 1? &e., ?j 19 &c. In other words, the value of 
jj>, q] is the same, whether it be obtained from the expressions for p, q in terms of 
the original variables, or by an analogous process from their expressions in terms 
of the new. 
Particular cases of (77-) are the relations 

[&,*<]=■— l, [&, £/]==[%, ^] = [4 S/]=0. ....... (78.) 

(See Theorems IV., V.) 

64. The relations (74.), (76.), (77-)> (78.) of the last article depend solely upon the 
form of the equations (73.), art. 62, which connect the new variables with the old ; 
and are independent of any supposition as to the equations which may determine 
either set of variables as functions of t. Let us now, however, introduce the suppo- 
sition that the original variables x l9 ... x n9 y l9 ...y n are determined as functions of t by 
the system of differential equations, 

l dZ l d/U /T \ 

*=a£ *=-s; (L) 

The relations just established enable us immediately to transform this system into 
another involving the new variables instead of the old ; for we have 

now '&4^ (see (74.), art. 63); 

dt uvji 

and if in the remaining term we substitute for x' p y\ their values from (I.), and for 
$i, fi their values % -p, it becomes 

dxj dyj dYii dr\i 

y (dZ dxj . dZ dyj 




dt\i dyj dtii 

which is equivalent to^-> if Z be supposed expressed in terms of the new variables. 
We have then ^—-j- + x~ 5 

df\i dr^i 

and, exactly in the same way, 

,__dV_dZ 

*' '<%. dil 

This result may be stated in the form of the following Theorem VIII*. If the system 

* This theorem, in its general form, is, to the best of my knowledge, new. But that case of it in which P 
does not contain t explicitly has already been proved in a different way by M. Desbovbs, who has* by means 

MDCCCLV. 2 U 
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of differential equations (I.) be transformed by the introduction of new variables 
n l9 ...n M connected with the original variables x l9 ... x n , y l9 ...y n > by the equa- 



»i? • • • ■»»? 



dP dP 

tions j-=^, -jy=ni, where P is any function of l u ... | w , y„ ...# n , which may also con- 
tain £ explicitly, then the transformed equations are 



(79 ) 



in which O is defined by the equation 

at 

and is to be expressed in terms of the new variables. (The substitution of the new 

dP 

variables in -^ is to be made after the differentiation. See art. 63.). 

Corollary. — If P do not contain t explicitly, -^=0 and <3>=Z ; so that in this case 

the transformation is effected merely by expressing Z in terms of the new variables. 

65. It follows from (77-), art 63, that tffg be any two integrals of the system (L), 
the value of [/, g] is the same whether it be derived from these integrals in their 
original form, or similarly obtained from the same integrals after transformation by 
the introduction of the new variables. And consequently if n integrals a 13 a 2 , ... a n of 

Tl (y\ mm 1 J 

the original system be given, which satisfy the-^-^ — - conditions [a*, aj=0, they will 

continue, after a normal transformation, to satisfy the analogous conditions, so that the 
method of finding the remaining integrals given in Theorem VII. art. 49, will also 
continue to be applicable. We had an instance of this in the case of the problem of 
central forces (art. 27.), where the above conditions were found to subsist after the 
transformation from rectangular to polar coordinates. (It will be shown presently 
that every transformation of coordinates is a normal transformation.) 

66. It was shown in Part I. (art. 18.), that if W be any function of x u x 2 , ... x n , 
x\, x 2 , ...x n (which may also contain t explicitly), the system of n differential equa- 
tions of the second order 



( 



™Y = *w 

dxi ) dxi * ^ 'J 

may be changed into a system of 2n equations of the first order of the form (I.), 

of it, deduced Jacobi's form of the method of the Variation of Elements (namely, the equations (68.), art. 54), 
from the similar form of Lagrange, in which the elements are the initial values of the variables. It will appear 
in the sequel that the extension to the case in which P may contain t, is of importance. If the expression were 
not already appropriated, I should have proposed definitively to call P the " modulus of transformation;" and I 
shall use this term provisionally in the present paper, not being able to suggest a tolerable substitute. After all, 
as the word " modulus " itself is used without confusion in very different senses according to the subject matter, 
there is, perhaps, no reason why a similar liberty should not be allowed in the use of the proposed expression. 
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art. 64, by putting y^-^r, and taking Z a function of x l9 &c., y 19 &e. and t 9 defined 

i 

by the equation Z==— W+S^a^), 

in which x 19 ...x n are to be expressed in terms of y 19 &c., x l9 &c. # . 

Conversely, a system of the form (I.) being given, it may be changed into a system 

of the form (80.) as follows: by means of the equations tf!=-p> let i/^ ...y n be ex- 
pressed in terms of x l9 &c, <r 15 &c. ; it follows from Theorem I. (art. 49.) that we shall 

dW 

have y^U • • •' (a.) 

£ 

dxi dxi x ; 

where W is a function of x 19 &c, *i, &c. defined by equation 

W=-Z+2 < (4y*), 
in which y„ ...#„ are to be expressed in terms of #i, &c, # 15 &c. The w equations 

#J= — ~= are then changed by (a.) and (b.) into the form (80.). 

"Qf30% 

On the Transformation of Coordinates. 

67. It has been seen in the preceding article, that we can always change the system 
of 2n equations of the first order of the form (I.), art. 64, into a system of n equa- 
tions of the second order of the form (80.). In this latter form the equations of 
dynamics naturally present themselves. 

Now in the case of the dynamical equations, x l9 x 2 ,...x n are the independent 
coordinates of the system (the word coordinates being taken in its most general sense), 
and when the equations are to be changed into the form (I.), the additional variables 

y l9 ...y n are defined by the equations -^T—y { . In this case a transformation of coordi- 

nates, in the most general sense, consists in taking n new variables % l9 g a , ...| w , con- 
nected with the original coordinates x 19 ...x n by n equations, which may also involve / 
explicitly. It is a well-known theorem, that the transformation of the equations (80.) 
is effected merely by expressing W in terms of the new coordinates l l9 ...% n9 and their 
differential coefficients £i, ... ?L, instead of the old ; so that the new equations are 

f dW\'_dW 

dg.) ~ dg t : * • ' * V 81 -) 

the proof of this theorem does not depend upon the form of the function W ; and we 
know also (see arts. 18 & 66.), that whatever be the form of W, these new equa- 
tions may be again transformed to a system of the form (L), by taking n additional 

* This theorem is a generalization of Sir W. R. Hamilton's transformation of the Dynamical Equations. 
See Part I. art. 18. 

2 u 2 
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variables n l9 ...jj n9 defined by the equations 

dW 



fli=- 



d? 



It is to be observed that this last transformation is not, in general, equivalent 
to expressing the original Z in terms of the new variables ; for the original Z 
is — W+2(%i) (art. 66.), and the analogous expression derived from (81.) is 
^W+2(g^), which is not, in general, equivalent to the former. It will be seen 
presently that the two expressions are equivalent when the equations connecting 
x l9 ... x n with \ l9 ... l n do not involve t explicitly.] 

68. Following the analogy of the dynamical equations, I shall adopt the following 
as the 

Definition of a transformation of coordinates. 

The original equations (I.), art. 64, having been changed into the form (80.), art. 66 9 
let | 13 g 3 , ... g n be n new variables connected with the n variables x 19 x 29 ... x n by n equa- 
tions, which may also involve t explicitly ; and let n l9 *? 2 , ...q n ben other new variables 

dW 
defined by the equations -^-=^ (where W has been expressed in terms of &, &c, £i, &c). 

By means of the 2n assumed relations, the 2n original variables x l9 ...x n , yi,-**y n can 
be expressed as functions of the 2n new variables &,...£,, n X9 ... n n - Let this substi- 
tution be called a "transformation of coordinates." 

It has been seen in the preceding article that the original equations (I.) are changed 
by a transformation of coordinates into a system of the same form, which however 
cannot in general be obtained by merely expressing Z in terms of the new variables. 
But we are not at liberty to assume (and it is not generally true) that a change of 
the system (I.) into another of the same form is a normal transformation (art. 62.). It 
has already been stated, however (end of art. 65.), that this is true in the present case ; 
a proposition which I proceed to establish. 

69. Every transformation of coordinates is a normal transformation. 

To prove this theorem, we have to show that every transformation of the kind 
described in the last article is also of the kind defined in art. 62; in other words, 
that it is possible to assign a function P, of l l9 ...£,, y 19 „.y n (with or without t), such 
that the given relations between x u &c, g w &c, which define the transformation of 
coordinates, shall be equivalent to the system of equations 

dP_ dP__ 

Take p =fe)yi+(^)y2+... + (^)y» ....... (82.) 

(the brackets indicating that x l9 ...x n are to be expressed in terms of £„•..£,; so that 
P is a function of f 15 ...L, y l9 ...y n9 with or without t according as the equations con- 
necting x 19 ... x H with f l5 ... l n do or do not contain t explicitly). Then P is the func- 
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dV 

tion required. For we have, at once, t~=^ ; also 



dW 
now the definition of ^ is ^=-^7-, where 

W ="~(Z)+% 1 +%2+...+^ 

(expressed in terms of f n ...§„, £1, ...§1 (art. 66.)) ; hence, putting Z (without brackets) 
for the original form of Z in terms of x l9 &cc. 9 y 19 &c., and observing that (Z) becomes 
a function of g' 15 &c, only through y„ &c, we have 

dW dZ dy x dZ dy^ dZ dy n 

d%~~~ dy x dg~dyz d%~~ '"'~~ dy n d£. 

J-*/ dXl JLii ^ 2 JL -Lit d0 °' n 

, dZ 
The two first lines of this expression vanish by virtue of the equations x'i=j- m ; and 

SlDCe ,^%),%)w l^Nv/ , . <%)*' 

X i— dt ^ //£. Si I" //£ 5iT-T 



A ^ #, * in € 2 ^-r — -r «$ 






■ Hub • fJ/Kobii) 

we have —?=-~^; 

so that we have, finally, 

dW d(x r ) d{x^) d(x„) 

which is evidently equivalent (see the expression (82.)) to 

dP 

the proposition in question is thus established, and may be enunciated as fol- 
lows : — 

70. Theorem IX. — Every transformation of coordinates is a normal transformation, 
of which the modulus * P is a function of% l9 ...% n9 y 19 ...y n (with or without t) given by 
the equation P=(^+(%+... + Wy. 

(the brackets indicating that x 1? ...x w are to be expressed in terms of% l9 ...% n ). (See 
arts. 62, 68.) 

This theorem will be made more intelligible by applying it to a very simple 
example. 

Let it be proposed, then, to transform from rectangular to polar coordinates the 
differential equations of any dynamical problem referring to the motion of a single 
material point whose mass is m. Let x, y, z be the rectangular coordinates of m 9 

* See note on Theorem VIII. art. 64. 
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1 

and u, v, w the variables conjugate to them ; so that, putting T=^m(x ,2 +y 2 +z ,2 ) 5 we 
have u=j 1 ,&e.; whence T=^(u 2 -fv 2 +w 2 ), and the equations of motion are x f =-^> 

u f =— -t-,&c, where Z=^(u 2 +v 2 +w 2 ) — U, and U is a given function of x, y, z, with 

or without t. 

Now let r, #, <p be polar coordinates of P, so that 

x=r sin cos <p, y =r sin sin <p, z=r cos #. 

Let u, v, w be the variables conjugate to r, t?, <p. Then the ordinary process of trans- 
formation would be as follows : — 

(1) to express x f , y', z' in terms of r, 0, <p, r', 0\ <p r 9 and thus transform T into a 

function of the latter quantities ; 

(2) to define u, v, w by the equations 

dT dT dT 

U =V V =W w =d?> 

and by means of these relations to express r', 0\ <p' in terms of w, v, w, r, 0, <p, so 
that x', y', z', and therefore, finally, T and Z, might be expressed as functions of 
the six new variables. 

Instead of this, let us adopt the method indicated by the theorem at the beginning 
of this article. 

We have then, for the modulus of transformation, 

P=(x)u+(y)v+(z)w, 

in which x, y, z are to be expressed in terms of r, 9 <p; so that the proper form of P is 

P=ur sin cos <p+ vr sin sin <p+wr cos 0, 

dV 
and the equations (corresponding to m=w. (art. 69.)) which define the new variables 

tM>,it>,are dV dJ > dp 

u =ip v =w w =w 

These give &=sin #(ucos<p+vsin<p)+wcos$, 

v=r cos 0(n cos <p+ v sin <p) —rw sin 0, 
w= —r sin 0{u sin <p— v cos <p), 

from which the values of u, v, w are easily obtained in terms of the six new variables. 
But in order to effect the transformation of T, we have only to square each side of 
these equations, and add them, after dividing the second by r 2 and the third by 
(rsin^) 2 ; we thus obtain 

V w 
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and the transformed equations of motion are therefore 

mr f =u, m(f=vr~*, mp f =w(r sin 0)~ 2 , 

d\J 
rnu'=r~ 3 (v 2 -\-w 2 (sm 6)~ 2 ) -\-m,-j-> 

mt/==(r sin 0) 6 w cos0+m-7p w/=-r-. 

The preceding process is not 5 of course, given for the sake of the result, which may 
easily be verified directly, but in order to illustrate the meaning of the theorem oh 
which it depends. It is hardly necessary to add, that if the problem involved the 
independent rectangular coordinates of any number of material points, the trans- 
formation to polar coordinates would be effected in the same way, by merely adding 
to P analogous terms for each point. 

71. Before proceeding further there is an important remark to be made. 

It has been hitherto assumed that the modulus of transformation P was a function 
of no other quantities than the 2n variables | 1? ... % n} y l9 ...y n and t. But if every step 
of the demonstrations of the theorems of transformation which have been given (art. 
62, &c.) be examined, it will be seen that they continue to hold good in the following 
more general form. 

1 ake ■^ == j v^i5 ?2? • • • bw? y\ 9 y^ • • • yw> Pi 9? ^\ • • •* v> 

where jo, #, r, ... are any functions of any or all the variables, old and new, with or 
without t. 

Let the equations connecting the old and new variables be, as before, 

with the condition that p, q, r, ... are exempt from differentiation informing these equa- 
tions. 

Then take T , = — (~irh with the following signification; (l)-^jr denotes the 

differential coefficient of P with respect to t, so far as t is contained explicitly, and 
also through the variables in p, y, r, ... ; that is to say, 

d(P) _dP dV t dP , • 

dt ~ ^+^r^ + ^ *?+••• 

(where p ! =j+j~~ x[-\-&c. &c, but this substitution is not to be made at this stage) ; 

(2) (~jt) denotes the result of substituting in the above expression the values of 

# 15 ...^in terms of % 19 ...%n 9 n l9 ...n n9 p 9 q, r,... from (83.), so far as — ^- contains y x ,... y n 
explicitly (L e. not involved in p, q, r, ...). 

Lastly, take 3>=(Z)+Y, 

where (Z) denotes the result of substituting in Z the values of the old variables as 
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given in terms of | 15 ... q 19 ... p 9 q 9 r, ... from (82.). We shall then have 



& =3T *="^' (84,) 

where O is in general a function of 

Sw ••• £» Hi, .- %., P, q, r, ... //, q 1 , r f 9 ... and # ; 

and the differentiations with respect to |„ % are performed only so far as those vari- 
ables appear explicitly in O. But after these differentiations , we may introduce the 
actual values of p 9 ... p' 9 ... in terms of the variables and their differential coefficients. 
It is obvious that the original variables will not, in general, have been eliminated 
from the system (84.) ; but of course the elimination may be afterwards completed*. 
Similar considerations apply to that particular case of transformation which we have 
called a transformation of coordinates (art. 68.). We have then 

and the relations connecting x x „.. x n with g 1? ... | rt , may contain p 9 p, r 9 ... ; so that 
(#j), &c. are functions of p 9 q 9 r 9 ... as well as of % 19 ... g n . 

We might have deduced the preceding conclusions from the following simple con- 
sideration. Since p, q 9 r, ... are actually functions of t 9 though unknown functions, 
we may imagine them to be known 9 and to be expressed explicitly in terms of t ; and 
then the case resolves itself into that of art. 62, &c., so far as the demonstration is 
concerned. But as a doubt might possibly have arisen whether any fallacy was 
involved in the circumstance that p 9 q 9 r 9 ... involve (when supposed to be expressed 
in terms of t) the arbitrary constants of the problem itself it seemed best to refer to 
the original reasoning; the most important part of which is that contained {mutatis 
mutandis) in art. 6. (For the " mutanda" see the beginning of art. 63.) It is then 
apparent that this circumstance is perfectly immaterial with reference to the con- 
clusions in question, though it may be important in other points of view. 

72. This being premised, we will proceed to an example of transformation more 
interesting than the former, namely, the 

Transformation from Jixed to moving axes of coordinates. 

Let x, y, z, u, v, w have the same signification as in art. 70, and let x 9 y 9 z. u 9 v 9 w 
be the new variables, where x 9 y 9 z are rectangular coordinates, referring to a system 
of moving axes of which the origin always coincides with that of the original fixed 
axes of x, y, z. 

Let the direction cosines of the new (moving) axes with respect to the old be 
\>, ^o> h ; K> pi, h \ K fa, h, thus-f : 

* The final equations in this case will not in general have the canonical form. 

f I do not know who first used this convenient way of indicating the nine direction cosines by a diagram, 
but I first saw it in one of M. Lame's works. 
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X 
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\ 


\ 


y 


jU'o 


^i 


fV 


z 


h 


h 


h 



where A , &c. are functions of t, which may be either given explicitly, or implicitly 
through the variables (see the last article). 

The modulus -of transformation P is found (art. 69.) by substituting for the vari- 
ables x, y, z in the expression xu+yv+zw, their values in terms of x 9 y, z ; we have 
therefore 

and then the three equations 

dV dV dV 

)U r Ty= V > Tz= W > £ lve 



u 



dx 
A u+^ v+vn 



u = x ^+ K Vm \" \ w 

v ^zXjU+^v+^w k whence I v ^f^+j^+j^w 
«i;=X 2 u+jM» 2 v+i'2wJ lw= V^+ "i^+ ^2^ 



? . 



(85.) 



Also we have (see Theorem VIII. art. 64.), since P is to be considered to contain t 
explicitly through X , &c., only, 

dV 

-^ — fox+ty +Kz)u + ((*<&+ (Aiy+faz) v+ {w+v$+v % %) w, 

in which expression the values of u, v, w in terms of the new variables are to be sub- 
stituted from (85.). Now if we put et> 09 'a> l9 co 2 for the angular velocities of the moving 
system of axes about the axes of x, y, z, respectively, so that 

it will be immediately seen that the usual relations between the nine direction cosines 
enable us to put the result of the substitution in the following form : 

/dV 



dt 



=a (yw~zv)^a 1 (zu~xw)-{-a) 2 (xv-~yu) 



The original differential equations 



, dZ , 

X= — -, U r 

d\x 



dZ 
dx 



? &c 



are then transformed into 



x = 



d$? f d<§ Q / , a A \ 

— , %'= — — > &c. (art. 64.) 

du dx 



where 

MDCCCLV. 



*=( z )-(f) 



2 x 
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Introducing the above value of (-7-), and omitting the brackets, we obtain for the 
system of transformed equations, 



x -=-~~-\-^ 2 y—oo l z, u = — -r — \-co 2 v— ca x w 

r &£k • f dZ 1 

y'=—-\-6j z — 6j 2 x, v = — — +&; z# — a; 2 & 

& = j — \-a) x x — o; ;?/, w; = — — ^-o){U — 6t> v 



(86.) 



dw ' x u «^ dz 

in which Z is supposed to be expressed in terms of the new variables. 

73. On the principles of the integration of this, and of transformed systems in 
general, I shall make some remarks hereafter. For the present, the following may 
be observed. If, in the transformation of the last article, we suppose the motion of 
the new axes given, then X , &c, and therefore also a , oo Y , a 2 , are given explicit func- 
tions of t. But if the motion of the new axes is only given by connecting it with the 
motion of the point m itself, then the above quantities are given functions of the 
variables and their differential coefficients. 

The most interesting case of the latter kind is that in which the motion of the new 
axes is assumed to satisfy the equations 



x y z y 



(«.) 



which express the condition that the instantaneous axis of rotation (of the moving 
axes) always coincides with the radius vector of the moving point m # . 

The radius vector traces, in fixed space, a certain conical surface. It also traces, 
with reference to the moving axes, another conical surface ; and we might always 
assume as one of the conditions defining their motion, that this latter should be any 
proposed surf ace ; that is, we might assume that the new coordinates x, y, z should 
always satisfy the equation <p(x,y, z)=0, <p representing any given homogeneous func-* 
tion. If to this last assumption we add the two conditions expressed by the formula 
(&>.), we further assume that the conical surface traced by the radius vector with refer- 
ence to the moving axes, rolls upon that traced infixed space. 

Suppose, for example, we assume for the equation <p(x,y 9 x)=0, simply z=0. This, 
with the conditions (a.), will express that the radius vector is always in the plane of xy, 
and that this plane rolls upon the conical surface traced by the radius vector in fixed 
space. We may then say that the plane of xy is the u plane of the orbit," and that 
the axes of xy, or any lines fixed with reference to them in their plane 5 are " fixed in 
the plane of the orbit-f\" 

* See Jacobi's first letter to Professor Hansen (Crelle's Journal, vol. xlii. p. 21). This letter appears to 
refer to some unpublished (?) results of Professor Hansen, which may possibly be similar to those of this 
article. 

t The student of elementary treatises is, I believe, always left to find out for himself what this expression 
means, or ought to mean. 
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Now since we are supposing the equations to belong to a dynamical problem, we 
have u=mx ! , v=my, w=mz', and if we substitute on each side of these equations 
the values derived from (85.), art. 72, we find easily 

relations which are true on all suppositions as to the motion of the axes ; but the 
assumptions (&.) reduce them to 

u=mx ! , v=my l , w—mz\ 

The further assumption %=0, which involves z f =Q, gives w=0, and also (by the 
equations (*>.)) ^ 2 =0. Thus the equations (86.) are reduced to 

du dx 

A r dZ i dZ 

U — ——-—!» V -. — — — - 

dv dy 


= -—■5 U= r-+#,W 0) Q V 9 

dw dz 

d*7 d7 

where in — and — , % and w are to be put =0. after the differentiation. It is to be 

dw dz l 

observed also, that the values of u,v 3 w above given reduce the three first of the 
equations (86.) to the form u=m —, v=m ^ w=m ~, from which it is evident that 

du dv dw 

Z=:±(u*+v*+w 2 )-XJ, 

where U is the original force- function, expressed in terms of the new variables. This 
depends only upon the conditions ('*.) ; but in the case now considered we have also 
v)=0. 

Let the origin of coordinates be the Sun ; m a planet disturbed by another planet 
m l whose original coordinates are x p y /? z y , and new coordinates x p y p z { ; also let 
x'+yM-z^r 2 , a*+y*+s?=r*, x *+y*+z*=r% tf+tf+z*=r% and (x,-x) 2 +(y-y) 2 
+ (z / «— z) 2 =S 2 ; we have evidently r 2 =r 2 , r 2 =r 2 , and 

(x,-x)»+( yi -y)»+(z,-z)»^ 
and xx / +yy / +zz / =^+^/ / +^^ / , 

We have then originally 



tt m, /I xx/4-yy. + zzA 



where ^=m+mass of Sun ; and it is evident that this expression preserves the same 
form when expressed in terms of the new coordinates, and also (which is essential to 

the validity of what follows), that ^, &e. are the same whether the differentiation 

Ji X Jt 
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be performed before or after the substitution for x ; y /? z/, in terras ofx p y l9 z r If then 
we put z=0 and 



R=m,{]-^>), 



I 



we shall have, from the fourth and fifth of the transformed differential equations 5 ^ 



r 3 dx 



7T1 

" r s dy 



> 



(87.) 



from which it is evident, that, assuming the motion of m { to be known relatively to the 
new axes, the variations of the four elements of the orbit of m which determine the 
dimensions of the orbit, its position relatively to lines Jixed in its own plane, and the 
time of perihelion passage, will be expressed in terms of the differential coefficients 
of R in the same way as if the plane of the orbit were fixed. But the motion of the 
node of the orbit upon the fixed plane of xy, and its inclination to that plane, must be 
determined by means of the last of the differential equations,, as follows : that equa- 
tion gives 

»,m— »o«>=(^)(s==0); 
or if we put — 12 for the term multiplied by mm { in the value of U given above, 



d f mu, , ^ \ 



and z is to be put =0 after the differentiation, which reduces the above to 

1 ° \dz 

Let e%+a%=o? 9 so that m is the angular velocity of the plane of the orbit about the 
radius vector; then (observing that u=mx l , &c.) we have 

30 ~~~~ V r m{xy ] —x ] y) 

whence v^ -^ua) l =~(xy f — x f y) 

and therefore m<%= 7=tT ==== w( jr )> 

which gives cc in terms of the four elements referred to above, and of t. And if we 
put / for the inclination of the orbit to the plane of xy, v for the longitude of the node 

* These equations (87.) have been obtained in a different way by Mr. Bronwin, Camb, Math. Journ. 
vol. iv, p. 245. 

f See below, art. 80. 
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referred to the axis of x, and (3 for the angle between the axis of x and the node, the 
usual formulae of rotation give 

t'=zoj cos 13-—^ sin j3 "| 

j/sin/=<a sin jS+a^cos/3 > (88.) 

fi ! z=a 2 — y'cos/ 

If in these expressions we put a> =- a, ^=^a, a> 2 =0, and call & the angle between 

the radius vector and the node, so that a/ cos|3— m x sin(3=rcos^ 5 a; sin^3+^icos^ 
=r sinS, we obtain finally 

i'=acos&, p f =a-i — > $= — a sin ^cot /. 

sin * ' 

in which the expression above given for a is to be substituted. 

The actual value of f -j~ ) i$mm j z l (-3 -""as h which (since z„ r p &c. are supposed 

given in terms of t) may be expressed in terms of the four elements first mentioned, 
and t. 

I propose to consider the transformation of the differential equations of the planetary 
theory in a more general manner in the following section. At present I shall add 
some remarks on normal transformations in general. 

74. Theorem. If 

/ dZt i dZi /on \ 

Xi =^; *="-s; •••••••••• (89 -) 

be a system of 2n simultaneous differential equations, where 

^ === J' v^D H^ ^2? 3/25 • • • P> $5 ^V • • • 5 *)» 

and p,q,r, ... are also explicit functions ofx 19 &c.,y 19 lkc. and t 9 but are exempt from differ - 
entiation in taking the differential coefficients _ , _ ; and if these equations be trans- 

dxi diji 

formed by a normal substitution of new variables | 15 &c, ^, &c. (art. 62, equation 
(73.)), then the transformed equations are, as in art. 64, 

£'=—4-—, '———-- — 

dm dm l dgi d%i 

in which Z is expressed in terms of | 19 &c, ^, &c, but the differentiations with respect 
to'ia ^ are performed before the substitution of these variables in jr>, q 9 r, &c.; in other 
words, p, 5, r, ... are still to be exempt from differentiation in forming the differential 
equations. 

This may be proved simply by repeating the reasoning of art. 64. The only dif- 
ference is, that in the term 

y \dZ d%j t dZ dyj 
J \dxj dv\i dyj dr^ 

the differential coefficients ^p> ^~ are now taken only so far as Z contains w p i/ f inde- 
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pendently oip,q, r, &c. ; and therefore the term represents the differential coefficient 
j~ taken so far as Z contains m independently of p, q, r, &c. The same reasoning 

applies to the corresponding term in the value of %. The theorem is thus established. 

It is evident that it may be combined with that given in art. 71j, where other func- 
tions analogous top, q, r, ... are introduced by the modulus of transformation P. 

If we call the form of the system of differential equations (89.) canonical when the 
differentiations of Z with respect to x x , &c, y l9 &c. are total, we might call it pseudo- 
canonical when Z contains functions of x 19 &c, y x , which are exempt from differentia- 
tion in forming the differential equations. 

In like manner, if we call a transformation of variables normal, when the differen- 
tiations of the modulus P (equations (73.),. art. 62.) with respect to f 2 , &c, y Y , &c. are 
total (as in art. 62.), we might call the transformation pseudo-normal when P contains 
functions of the variables which are exempt from differentiation in forming the equa- 
tions of transformation (as in art. 71-)- 

Adopting these designations, we may enunciate the following general theorem of 
transformation : — 

Theorem X. — If a pseudo-canonical system be transformed by a normal or pseudo- 
normal substitution, the transformed equations are also pseudo-canonical, and may be 
formed by the rules applying to normal transformations of canonical systems, provided 
that the functions which are originally exempt from differentiation with respect to 
the variables, be continued exempt to the end of the process ; but if such functions 
occur in the modulus of transformation P, they are subject to total differentiation 

d>P 

with respect to t in forming the term -rr* (See art. 71 •) 

[With respect to this theorem there is one important remark to be made. If w, v 
be any two functions of x Y , &jQ.,y x , &c. (with or without p, q, r, ... and t), the equation 

'du dv du dv\ /du dv du dv\ ( 

dtjidxi dxi dyj \drn dfa d& dfa) ^ * •' 

is now only true on condition that the substitution of the actual values of p, q, r, ... 
in terms of the variables be not performed till after all the differentiations.] 

75. The theory of the variation of elements affords an interesting example of the 
theorem given in the last article. Consider the following system of differential 
equations, 

Xi ~ d Vi + dyS y i — w ~dx~dx i 9 • • • ■ ( 90 «) 

dZ dZ 

where in ^ j- the differentiations are total, but H is supposed to contain functions 

of x x , &c, y x , &c, which are exempt from differentiation in forming the above equa- 
tions. The system 

,__dZ ,_ dZ 
1 dyi y* ~dxi 
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is canonical. Let us assume then that a complete set of normal integrals a l9 ... a n9 
b 19 ... b n of this latter system is known, so that we have 

tf«=ft(#i, &c. 5 y 19 &c, *)> h=-xl^ &c, #„ &c, f). 

The assumption of these last equations to represent the solution of the complete 
system (90.) is simply a transformation of variables, (a 19 &c 5 b l9 &c. being the new 
variables) ; it is also a normal transformation, since the equations connecting the 
new and old variables may be put in the form (see Theorem VII. art. 49, and art. 62.) 

dX_ d JL-. h 

where X (the modulus of transformation) is a function of x l% ... x M a l9 ... a n9 t. The 

function T" of art. 62 is now obtained by expressing — -^ in terms of a 15 ... & w ... ; 

dX. 
but since Z is — - -rr expressed in terms of x l9 &c> 9 y 19 &c., it follows that when Z is 

expressed in terms of the new variables a l9 &c., b i9 &c, it becomes identical with "¥". 
Now if the process of art. 63 be followed, mutatis mutandis, it will be seen that in 
the present case we obtain 

in which expression the first two terms destroy one another, and the remaining term 
is evidently the differential coefficient of Q with respect to b i9 taken so far as Q con- 
tains b t independently of those functions which were exempt from differentiation in 
forming the original differential equations (90.). Similar reasoning applies to the 
expression for b' { . 

As this result will be useful, I shall enunciate it separately as 
Theorem XI. — If the original system of differential equations be formed by treating 
certain functions, p 9 q 9 r, ..., contained in the disturbing function Q, as exempt from 
differentiation with respect to x 19 !kc 9 y l9 &c, the equations which determine the 
variations of any set of normal elements a l9 &c, b l9 &c. 

, d£l v dVt 

on condition that p, q, r, ... be treated, in forming these equations, as exempt from 
differentiation with respect to a 19 &c, b 19 &c. 

[It is important to recollect, that after these equations are formed, p 9 q 9 r, &c. are 
to be expressed in terms of a 19 &c, b l9 &c, and in the integration of the system 
a 19 &c, b 19 &c. are to be treated indiscriminately as variables, whether they originally 
entered through p 9 q 9 r, ... or not]. 

The Theorem XI. may also be immediately obtained from the general equations 
(E.) of art. 52 (in which it is to be remembered that Z includes the disturbing 
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function). The above method of deducing it is given as an additional illustration of 
the general theory of transformation. 

If, instead of the normal elements a lt &e., b 19 &c, we employ any other elements, 
c 15 c 2 , &c.j which can be expressed as functions of the former, the formula (69.) of 
art. 54 will still be applicable, with the condition that p, q, r, &c. are exempt from 
differentiation with respect to c l9 c 2 , &c. 

Section VII. — On the Differential Equations of the Planetary Theory. 

76. The differential equations which determine the motions of a system of mutu- 
ally attracting material points relatively to one of them, do not, as is well known, 
naturally present themselves under the canonical form (I.), art. 49. It is possible 
indeed to reduce them, by different artifices, to that form; but it seems doubtful 
whether any practical advantage is gained by doing so. When the ordinary method 
is followed in the case of a planetary system referred to the sun, there is a distinct 
disturbing function for each planet ; but it is easily seen that the usual expressions 
for the variations of the elements hold good, not merely for each planet on the hypo- 
thesis that the motions of the rest are known, but as a complete and rigorous set of 
simultaneous differential equations involving all the elements of all the orbits, and 
their differential coefficients with respect to t (and containing of course also t ex- 
plicitly). It does not appear that we are practically farther from the attainment of 
the rigorous integration of this system, than we should be if it had the canonical form, 
as it might be made to have if it were derived from an original system of that form ; 
and so far as the development of the disturbing functions is concerned, the most 
troublesome part of them, which is that depending on the mutual distances of the 
planets, is not likely to be got rid of by any conceivable artifice. 

However this may be, all that I propose to do in the present section is to take the 
original differential equations in their ordinary form referred to rectangular axes 
passing through the sun and parallel to fixed directions, and then to exhibit in a 
general manner the effect of a transformation to new rectangular axes, still passing 
through the sun, but changing their directions in space according to any arbitrary 
law. 

77* Let M be the mass of the sun, and m, m p m jP &c, the masses of the planets ; 
and put M+m=^, M+m^^, M+w w =^ ... M+m (i) =:{jt, (i) . And, referred to the 
original axes, let x,.y, z be the coordinates of m, x /5 y p z. of m p ..., x m y (0J z m of m (i) . 
Also let R, R p ... R (0 have their usual significations, so that 

= 2 ) m d) aft ffl(xx ( <) + Jj (i) + zz m ) 



the summation extending to all the planets except m. 

Then if we put xf +yf +zf i) =rj )f the original differential equations of the second 

order are such as x" i) +^=^. Let u w , v (0> w (i) be the variables conjugate to 



■Q) tt-X-(t') 
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x (t)> Ym %• Then instead of the original system of differential equations of the second 
order, we have the following system of the first order : 

jrr jr? jrw 

^"^ y «"~^ 5 Z(0 ~3^ 

in which Z^ss^u^+vfo+wf^-m^^+R^). 

These equations are not of the canonical form, because R (i) is not the same for all the 
planets. But it is easy to put them in the pseudo-canonical form (art. 74.)? a process 
which is not necessary, but saves trouble by bringing them under the operation of the 
general rules of transformation established in former articles. 
In fact, if we take 

T r 



\ 2m (i) r«) / 



) x (y)+yfflyo-)+z (0 z (i) x^x^+ycoy^ + z^z^ 



\ 1 C/) 



■HA 

*(0 



_ 1 \ 

^((X(«~x . ) )2 + (y (^ • ) ^y (y) ) 2 +( Z(i) ^z (y) ) 2 )^ , 

where the summation in the first term extends to all the planets, and in the second 
to all their binary combinations, and a horizontal line placed over any letter indicates 
exemption from differentiation, we shall have 

v dZ i d ®> w — . dZ dQ > (91 ) 

with similar equations for y (i) , z' w , v' (0 , w (i) . 

[The terms ^jj7 &c. are only written for the sake of uniformity, being really =0, 

since Q does not involve u (0 , v (0 , w ({) .] 

In these equations Z and Q are the same for the whole system, and the differentia- 
tions of Z are total; but those of Q are restricted to the quantities not marked by the 

horizontal line, so that ^~- is really the same thing as — ^(o^T* 

78. Let us now refer the whole system to new (moving) rectangular axes, whose 
position at any instant with respect to the original axes is defined, as in art. 72, by 
the variable direction-cosines X , &c. Let x {i) , y {i) , % {i) be the new coordinates, and 
%)> %)> w (i) the new variables conjugate to them. The transformation will be effected, 
as in art. 72, by taking for the modulus 

P=2((^ (0 +X I y (0 +X a *»)u W) + (fMta+l*iy<o+l^(o) v (o+ ( Vto+'i^o+^o) w (i) ), 
and the result will be as follows : put 
n=Q+5 2(* (0 v w — y w ^ 

MDCCCLV. 2 Y 
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(where Q is expressed in terms of the new variables, and &> , w l9 a 2 (the angular velo- 
cities of the moving system of axes about the three moving axes themselves) are 
marked with the horizontal line to show that these quantities are exempt from diffe- 
rentiation in forming the following system of differential equations, though they may 
be functions of the variables) ; then the system (91.) is transformed into 

, __ dZ d£l , dZ d£l . 

w du(i)*dua) ' (i) ~"~~ dotty dxtff ' ' ' ' * * * ^ *' . 

with similar equations for y u &c. 

In these equations Z is to be expressed in terms of the new variables ; and it is 
evident from the original form of Z and Q, that when so expressed, these two quan- 
tities are the same functions of the new variables that they were of the old, and in- 
volve (see art. 74.) the quantities exempt from differentiation in the same way*. 
Thus the transformed system (92.) contains no terms explicitly depending upon the 
motion of the axes, except those introduced by the three terms multiplied by w os a> 19 &> 2 
in the value of CI given above ; and the addition of these terms constitutes the only 
difference between the form of the old and of the new system. 

79. We may now apply the method of the variation of elements to the system (92.) 
as follows : — 

The system obtained by omitting the disturbing function O, namely, 

, dii , d/i t dju 



, dZ , dZ , dZ 



v.. .... . (93.) 



is canonical, and consists simply of the aggregate of the equations representing, for 
each planet, undisturbed elliptic motion about the sun ~f (relatively to the new axes 
of coordinates). 

The integrals of these equations may therefore be expressed in any of the usual 
forms. We will suppose that the elements chosen are 

a, e, w, (g), /, v, 

with significations corresponding to those given to the same symbols in art. 55. 
These letters unaccented will apply to the planet m, and a p e }3 &c, a in e in &c, a (i)9 e (i) , 
&c, to the planets m t , m u , wi (i) , &c. 

The definitions of the elements a, e, w, &e. are their expressions in terms of the six 

* Since the direction cosines \ , &c. are exempt from differentiation in forming the equations connecting the 
old and new variables from the modulus P, they continue exempt throughout. (Theorem X, art. 74,) Hence 

we ave x */+yy/+ z ^=(V+^+V^^ 

and similarly for the rest. 

fu % + v% -j- w 9, (&m\ 
t Z=S q- J, the summation extending to all the planets. 
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variables x 9 y 9 % 9 u 9 v 9 w* 9 and t 9 and the same expressions continue to be their defi- 
nitions when they become the variable elements of the disturbed motion. 

Now the general formulae for the variations of the set of elements here chosen 
have already been given in art. 55 ; for it is evident that the process in the present 
case would merely be a repetition, for each planet, of the process there employed. 
Here however we are to use in every case the disturbing function O given in art. 78 ; 
but if we observe the effect of the marks of exemption, it will be evident that, for the 
planet m 9 the only effective part of O is 

and similarly the effective part of Cl for any other planet will be given by suffixing 
the corresponding number of accents to x 9 y 9 z, u 9 v 9 w, m 9 R. 

Now if we put A , A 15 A 2 for the terms multiplied respectively by w 09 a l9 a> 2 in the 
above expression, we have, by the definitions of the elements, 

(A^+A?+A 2 2 ) i =m x /^a(l-e 2 ) 

A 2 = — 'm\/^a\l — e 2 ) . cos /, A x =ms/pa(\ — e 2 ) . cos v sin / 
A = —m\/pa{ 1 —e 2 ) . sin v sin /, 
so that the disturbing function, so far as m is concerned, becomes 

—m{R-\-\/(Aa(l — e 2 ). (oi sin vsm/—^ cos^sinz+^cos*)}. .... (94.) 

Consequently, since the expressions in art. 55 were obtained by taking -mR for the 
disturbing function, we have merely to add to them the additional terms derived 
from the part added to R in (94.). Performing the differentiations, and omitting 
afterwards the symbols of exemption over ^ *>i> ®2 which cease to be of any use, we 

obtain, after obvious reductions, the following simple results: if -^> ~^> &c. represent 

those parts of the differential coefficients of a 9 e 9 &c, with respect to t 9 which depend 
upon the motion of the axes of coordinates, then 



dt dt 



~fi = _. = tan ~(m\ cos v— &) sin v) —<*>$ 

at at ^ 



dt 



(»„ cos y+^i si n "J 



d 



V 



(95.) 



= — <y 2 + cot / (^oSinv — QJ l cosv) 
where it is evident that we may write s instead of (g) (see art. 55.). 

* Not their expressions in terms of w, y, z, oo\ y\ z'; for though these are equivalent in the undisturbed equa- 
tions, they are not here equivalent in the disturbed equations, and therefore the general theory, which assumes 
the former mode of expression, is not here applicable to the latter. 

2 Y 2 
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The complete variations of the elements are then found by adding to the terms 
just written the expressions given in art. 25. 

It is easily seen that the expressions (95.) might have been deduced from geome- 
trical considerations alone, if we had been at liberty to assume beforehand that the 
mechanical and geometrical parts of the variations might be calculated separately; 
the former as if the axes were at rest, and the latter as if there were no disturbing 
forces. It would not, I believe, be difficult to establish by a priori and simple 
reasoning the validity of such an assumption, and then the above results would only 
serve as a verification of the method which has been employed to obtain them. 

80. In order however that no obscurity may rest upon the interpretation of the 
formulae obtained in the last articles, it is necessary to consider the physical (or rather 
geometrical) meaning of the elements a, e, &c, which we have so far only defined by 
means of their expressions in terms of the variables x, y, z, u, &c, and to ascertain 
what relation they bear to the elements similarly defined by means of the original 
variables x, y, &c, which refer to axes whose directions are invariable. 

The relations between the variables ((85.), art. 72.) give immediately 

# 2 +# 2 + # 2 =x 2 +y 2 + z 2 

w 2 -]_ tf -J- w 2 — u 2 + v 2 + w 2 

ux-\-vy-$-zw=iix-\-vy-\-zw ; 
and if we put yw—zv~A, zu—xw=B, xv—yu=C 

yw— -zv=A, zu— xw=B, xv— yu=C, 
we find, by virtue of the relations ^ a — ^6o 2 =X , &c, the following equations : 

A=\ A+[* B+i> C 

C=k 2 A+p 2 B+p 2 C 
and ^ 2 +5 2 +C /2 ==A 2 +B 2 +C 2 . 

Now A(=yw— zv=m(yz'— zy')) is the projection on the plane of yz of the areal 
velocity of m (relative to fixed space) multiplied by the mass, and B, C have analogous 
meanings; hence it is evident from the above equations that A, B, Care the pro- 
jections on the three moving coordinate planes of yz, zx 9 xy of the absolute areal 
velocity of m relative to fixed space, multiplied by the mass. (The projections of the 
areal velocity relative to the moving axes would be yz ! —zy', &c, which are not pro- 
portional to yw— zv, &c, since u, v, w are not the same as mx ! , my 1 , mz f , except on a 
particular hypothesis as to the motion of the axes. See art. 73.). 

Inasmuch as the definitions of the elements a, e, i involve the variables only in the 
forms x 2 +y 2 -\-z 2 , u 2 -\-v 2 -\-w 2 , A, B, C, it follows that these three elements are 
respectively the semiaxes, excentricity , and inclination to the plane of xy, of the absolute 
osculating ellipse of the orbit in fixed space. Thus the instantaneous ellipse, relatively 
to the moving axes, is of the same dimensions and in the same plane as the true 
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osculating ellipse ; and it only remains to show that it coincides with the latter in 
position, for which purpose we must prove that it touches the true orbit. (It does 
not in general touch the relative orbit, because mx\ wy\ mz 1 are not in general the 
same functions of the elements and t that u 9 v, w are.) 

81. If we suppose the coordinates x 9 y 9 % of m expressed in terms of the elements 

7 

a, e, &c. and t 9 and denote by -^> &c. their differential coefficients taken so far as t 

fi ¥ fin ti ? 

appears explicitly in these expressions, then — 9 -j& -n are proportional to the direc- 
tion cosines, relatively to the axes of x, y, z 9 of the tangent to the (relative) instan- 
taneous ellipse. And therefore the direction cosines of this same tangent referred to 
the fixed axes of x, y, z, are proportional respectively to 

. dx , - dy , ~ dz dx , dy . dz dx , dy . dz 

^Tt+Ht+^J? ^Tt+^i+P'dt' v >Tt +v 4t+ v >1t 

On the other hand, the direction cosines of the tangent to the absolute instantaneous 

ellipse, referred to the fixed axes, are proportional to -J> X J (the differential 

coefficients of x, y, z, taken so far as t appears explicitly in the expressions for those 
variables in terms of the elements of the absolute ellipse). The identity of the two 
tangents will therefore be established, if we can show that 

dt lt' Kl Tt +X2 It 96cQt 

dx 

Now -j t I s that part of x f which does not depend upon the disturbing function ; i. e. 

(equations (92.), art. 78.) we have identically 

dx dZ dy__dZ dz _dZ-. 






dt du dt dv dt dw 

din dx dZ dv dZ dz dZ 

, in like manner. = , 4-=-t-> -r=-r L > 

dt d\x dt dv dt dw 

where Z is the same as before, but expressed in terms of x, y, z, u, v, w, instead of 
x.y y %)U 9 v 9 w. Now let the latter set of variables be expressed in terms of the former 
by the formulae of art. 72, and we have 

dZ dZ du dZ dv dZ dw 

du du du l dv du. l dw du "> 

du dv dw , . ^^ x 

but SI—*" du~ =x » du~ =x * (( 85 -)' art 72 -)> &c -> 

!• r dx ^ dx . dy t dz a 

and therefore lft =Xo ^+ x ^+ x vp &c *' 

as was to be proved. 

82. It follows, then, that the mode of treating the problem adopted in the pre- 
ceding articles is equivalent to representing the motion of each planet by means of 
the true osculating ellipse of its actual orbit (relatively to the sun) in fixed space. 
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The definitions of all the elements (relative to the moving axes) in terms of the six 
new variables x 9 y, z 9 u, v 3 w, have the same form as those of the corresponding 
elements (relative to the fixed axes) in terms of x, y, z, u, v, w. The two relative 
elements a, e are the same as the corresponding absolute elements ; / is the inclina- 
tion of the plane of the ellipse to the moving plane of xy, and v the longitude of the 
node reckoned from the axis of x; and since the place of the body in the ellipse is 
evidently the same, the relations between the remaining elements wand (e) (org) and 
the corresponding absolute elements are purely geometrical. 

Comparing these results with those of art. 79, we see that the independence of the 
formulae for the mechanical and geometrical variations of the elements of the true 
osculating ellipse is completely established. 

83. In all that precedes, the three variables # , a> 19 qj 2 (the angular velocities of the 
system of moving axes about the axes themselves) are entirely arbitrary ; they may 
be either explicit functions of t, involving only determinate constants, or they may 
depend in any way upon the relative or absolute elements of the orbits of any or all 
the planets, and their differential coefficients with respect to t. In the case in which 
the expressions for &/ , a> 19 & 2 involve only the relative elements, when these expressions 
are introduced in the formulae (95.), art. 79, and these formulae completed by the 
addition of the terms in art. 55, and when the corresponding sets of equations are 
formed for each planet, we obtain a set of simultaneous differential equations in- 
volving all the elements of all the orbits, and their differential coefficients with 
respect to t. The integration of these equations would determine all the elements 
as functions of f,' and thus the motion of all the planets, relatively to the axes of 
coordinates, would be known. Lastly, the motion of the whole system, relatively to 
fixed space, would be found by integrating the system of equations 

U=co cos X— oj 1 sin X 
N'sinL^sinX+^cosXl • • ♦• • ■• • • (96.) 



X f ~^ 9 — N'cosL 



'2 



■ 



where &> 09 &„ & 2 are now given functions of t, and L is the inclination of the plane of 
xy to that of xy, N is the longitude of the ascending # node of the plane of xy, 
reckoned from the axis of x, and X is the longitude of the axis of x, reckoned upon 
the plane of xy, from the node, in the direction of positive rotation. 

In the case in which <y , &> u &> 2 cannot be expressed in terms of the relative elements, 
the integrations which determine the relative motion of the system cannot be sepa- 
rated from those which determine the position of the axes in fixed space ; but the 
equations (96.) must be considered simultaneously with the other differential equa- 
tions of the problem. 

* Ascending relatively to a positive rotation, i. e. from $ to y. 
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Application to the Problem of Three Bodies. 

84. I propose to exemplify the preceding principles by applying them to the trans- 
formation of the equations which determine the motion of three mutually attracting 
bodies, considered as material points. Let them be called the sun and two planets, 
and let the origin of coordinates be placed at the sun, and the notation be the same 
as in arts. 77-79, so that M, m, m l} are the three masses, M+m=^, M+^=^ &c. 
Let the elements be chosen as in art 55, and longitudes be measured, as before, 
along the plane of xy (from the axis of x) as far as the node, and then along the plane 
of the orbit. 

Putting for convenience mR=0, m^^flp we have 

n==mm / {(r 2 +^--2rr y cosx)""^--"^r 2cos %} 
rL l =mm i {(r 2 -\-r?--2rr l CQ$ %)~*— /;r~ 2 cos %}, 

where % is the angle between the two radii vectores. Let I be the mutual inclina- 
tion of the planes of the two orbits, and let the angular distances, on the planes of 
the orbits, between their ascending nodes on the plane of xy and their line of inter- 
section, be respectively v 9 v t \ so that 6—v—v 9 d j —v l —v l are the angular distances of 
the two radii vectores from the line of intersection (which we may call, simply, the 
line of nodes). We shall then have 

cos 3£= cos (0—y—y) 008(0,— u^— v / ) + sin(5— y— • j>)sin(^«— y,— ^)cosI"| ^ . 

cos I = cos / cos /,+ sin / sin /,eos (v / — ») # J 

and y, v t are functions of /, / /5 v,— v, determined by the equations 

cot v sin (v / — v) = — cot / y sin i+ cos (^— • v) cos f 

T * * * • » 

cot w^sin^— ») = cot / sin /,+ cos (^— *>)cos / ; . 

Now considering H as expressed, on the one hand, in terms of r, r n 0, p v, v n i 9 $ p and 
on the other in terms of all the elements and t, we have rigorously, as may easily be 
proved in the usual way-f~, 

<$ afs ' d^ dij ds t "* dv*l 



(98.) 



* The arrangement referred to here and in the following 
articles will be made clear by the accompanying diagram. 

f Since the values of r and in terms of the elements and 
fare of the forms 

we have .^.j. ___=(), T -f--=l; 

as dm as dm 

from which the equations in the text follow immediately. It 
may be as well to remind the reader who may happen to 
recollect the note to the Astronomer Royal's tract on the 
Planetary Theory (p. 91, ed. 1831), that that note refers to 
a different way of measuring longitudes. 
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with similar equations for Cl r In what follows, I shall, as an abridgment, employ 
the symbol E to denote the operation j+-t-> and in like manner I shall put E, for 

£+4-> so that 

Efi=$+f , &c. . ■ (99.) 

as d'UT v ' 

Since r and d> when expressed in terms of the elements and t, do not contain /, /„ v, v n 
we have 

d£l dH dco&x d£l d£l dcos% 

dv d cos % dv di d cos % d\ 

and dcosx _.__( l idv\ dcos X _dv i dcosx 

dv \ ~ dv) di dv d6 t 

+ sin/ sin / y sin (0 — y— v) sin (0,— a,— v) sin (v,— v), 

in which expression the values of —-> -^ are to be obtained by differentiating the equa- 
tions (98.). In like manner we find 

c?cos% dv dco&x dv t dco&x 

di di di di *di { 

+ sin (0— u— v) sin (0 / — y y — *,)( — sin/ cos / y + cos / sin /, cos (vy—")). 
Analogous expressions may be found for — J> —J- 

UV i Ult 

85. Hitherto we have assumed nothing concerning the motion of the axes of 
coordinates. Let us now however take as a first assumption that the plane of xy shall 
always pass through the line of nodes. This implies the conditions # 

and consequently 1=^— /. 

Introducing these conditions in the expressions at the end of the last article, and 

* The legitimacy of the following processes will be apparent to the reader who shall have followed the 
general reasoning of former articles, though probably not to others. In either case it may be useful here 
briefly to recapitulate the principles now to be applied. The results of art. 79, and in particular the expressions 
(94.), (95.), were established independently of any assumption as to the values of w , w v w 2> which are perfectly 
arbitrary. We are therefore at liberty to assume that w , w p cv 2 are such functions of t that any three func- 
tions of the variables shall constantly =0. Thus the first assumption made in the text is that v f — y shall con- 
stantly =0. But since such assumptions may cause (as in this case) certain elements to disappear from the 
expression of £1, it is necessary to perform the differentiations of Q, with respect to such elements first ; the dif- 
ferential coefficients may or may not vanish, on afterwards introducing the assumptions ; if not, we find ex- 
pressions for them in terms of differential coefficients with respect to other elements which do not disappear ; so 
that instead of differentiations which ought to be performed before the assumed conditions are introduced, we 
have finally only such as may be performed afterwards. The expression " disappear" does not, it must be 
observed, necessarily mean the same as "vanish." Thus the condition Vf— y=0, causes v and v t both to dis- 
appear from O, but does not of course imply that they both vanish. 
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observing that -^-: 



EQ, &c. (see equation (99.)), we have evidently 

d£l 7 ' , /dv\ \ /dv, 



dv ' 

m 

di 



-(•+®)«°-@>p 



i)™-(%y,°-fv 



(7 \ 1 

j\, &c, namely, the values of ^>.&c., 

which correspond to the assumed conditions ^— y=0, t>=t> y =0. 

Now in any spherical triangle of which a, b, c are the sides, and A, B, C the oppo- 
site angles, if a be considered as a function of c, A, B, by virtue of the equation 

cot a sin c= cot A sin B+ cos c cos B, 
we have by differentiation 

dec -w~%\ 

-r= (sin a) 2 (cot a cot c+ cos B) 
= cos a sin a cot c+(J. — (cos a) 2 ) cos B= cos B+ cos a cot C sin B 

da / sin « \ 2 sin B sin « sin B 



dk. 

da 
dB 



/sinflV 
u — \sin A/ 



sin c 



sinC 



( sm a > cot A cos B+ cot c(sin af sin B= sin a cot C ; 



smc 



and if the sides of the triangle be indefinitely diminished, these expressions become, 

in the limit, 

da sin(B + C) da ~ da^ ^ 

~dc sinC 9 dK~ 9 dB~~ 9 

provided the angle do not vanish. 

If these results be applied to the triangle of which the sides are v, v p ^ — v, and the 

jY &c. are as follows : 



opposite angles or— i l9 /, I, it is easily seen that the values of 

shii, 



dv 
dv 



sm* 



I. 



sinT 



dv 



sinT 



dv { 
dv 



__ sin * / dv\ 
sin I \dvj 



sin* 
sin I 



dv 

dt 



dv^ 
di 



dv 



-' =o, 



dv _ m 

a) ~ V*7 ~ W 

provided I be not =0. We have therefore 



d£l 
dv 



EQ + ^i(sin/ / .EQ+sin/.E ) Q) 



da 



dv d\ 



(100.) 



exactly in the same way we find 



MDCCCLV. 



■%}= — E,Q,— sn(sin i .E,Q,+ sin <,. EQ,) 

da i da l 

di l mmm " dl . 

2 z 



e • 



(101.) 



(102. 
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The variations of the elements will now be found by introducing the above values of 

d£l dQ,, dQ, d£l. 

■ . 5 1, , L 

dv dv t di d\ l 

in the expressions given in art. 55, and completing these expressions by the addition 
of the terms (95.), art. 79. But the angular velocities &> , v 19 # 2 are no longer wholly 
arbitrary, since we have made one assumption concerning the motion of the axes, 
which implies one relation between these quantities and the elements. In order to 
determine &> , & l9 & 2 completely, it will be necessary to make two more assumptions ; 
but first we will investigate the relation already implied. 

86. The complete expressions for v\ v fl obtained in the way mentioned in the last 
article, from arts. 55 and 79, may be written in the following form : put for brevity 

and put ft*a-S for rc, and ftfaf-% for n t ; then 

/== — = — -i — h cot /.(a> sin v— & x cos v) — % 

i JJC% 

v= — -= — t j + cot iJ&q sin v,— «! cos v.)— w 2 

* p { sin *, di t * ' v ° ' * " 2 , 

and if the latter equation be subtracted from the former, and the conditions 
be introduced, the result is easily found to be 

, v . T sin* ^O, , sim, ^Kl 

(» smy— ^cosvjsmls— ■j^+-y i -^j. .... . (103.) 

This is the relation between a> , ^ and the elements and /, implied by the one 
assumed condition that the plane of xy passes through the line of nodes. The angular 
velocity &> 2 of the system of axes about the axis of z, is so far left, as it evidently 
ought to be, perfectly arbitrary, 

Q and Q / are now functions of the following elements # only : — 

a > e> 00> *r, % e n ( B i)> w v *> "• 
And we now have 

cos %= cos (d—v) cos {6 t — v) + sin (6—v) sin {6—v) cos I. 

87. The complete expression for / f , derived from arts. 55 and 79, is easily put in 
the form 

/ ' = ^^in^{"^ + ( 1 "~ C0S/ ) EQ — (»oCOSJ>+»iSin*) . 

and, on introducing the assumptions of the preceding articles, this will be found to 
become, after simple reductions, 

p sin I./= — (cos I.EQ+E^) — (iy cos y+^sin vjp sin I ; 

* On the difference between (s) and s see above, art. 55. We cannot strictly call a a function of a and s. 
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similarly, we find 

^ / sinI./=cosI.E / Q / +EQ / — (^oCosy+^sin^^sinl ; 
and since 1'=/— •/, we obtain from these 

p^sinU^^^osI.E^Q^EQ^+^^cosI.EQ + E^Q), . . . (104.) 
an equation which may be transformed as follows : — since 

i=°> |=0 (art. 79), 

n dSl dH 

we find from art. 55, y=-^-+^— = EQ, 

and similarly p t =:Eft l9 whence it is easily seen that the above equation may be written 

in the form 

(^cosI)'+ j pEQ i +^E i Q = ; . (105.) 

If we investigate, from the above expressions for / and /, the values of (p sin/)' and 
(ftsini^wefind 

sinI.(josin/) f =--(cos/ r E + cos/.E / )Q--j9Cos/sinI(^ cost'+A; 1 sinp), 
sinI.(p / sin/ / ) / =(cos/.E / +cos/ i .E)Q i — ^cos^sin I^cos^+^sin^), 

and, adding these equations, 

sin I.(psin/+P/Sin/ / )'=(cos/ / .E+cos/.E / )(Q / -— Q) 

— (p cos/ +JP, cos / y ) sin I . (» cosy +^ sin v). . . . (106.) 

With respect to this equation and (103.), it may be observed that ^cosv+^sin v 
is evidently the angular velocity of the plane of xy about the line of nodes, and 
^cosv— Go smv is the angular velocity of the same plane about a line in itself per- 
pendicular to the line of nodes ; or;, which comes to the same thing, the angular 
velocity of the line of nodes itself in fixed space, estimated perpendicularly to the plane 
ofxy. 

88. The position of the plane of xy at any instant has been so far left arbitrary, ex- 
cept that it has been subjected to the condition of passing through the line of nodes. 

As a further assumption, that which most naturally presents itself is, that the plane 
ofxy should always coincide with the principal plane. By the principal plane I mean, 
of course, that on which at any instant the sum of the projections of the areal velo- 
cities (multiplied by the masses) of the two planets about the sun, is a maximum ; it 
evidently always passes through the line of nodes, and would be the invariable plane 
if the disturbing functions vanished. 

To determine the position of this plane we have (see art. 80.) to express the con- 
dition that / and /, are always so taken, subject to the equation /,-— / = I, that the ex~ 

pression m\/jjja{\ — e 2 ).cosi-\-m i \/(jb i a i (l—ej i ).cosi i shall be a maximum. We will 
put a for the value of this expression, so that using p and p i with the same meaning 
as before (see art. 86.), we have 

c= i pcos/+p / cos/ / , 
2 z 2 
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and the required condition of a maximum will evidently be 

^sini+j^shw^O'*; . . . . . . . ' . (107.) 

adding the squares of these expressions, we obtain 

G*=f +^ + 2^008 1, . (108.) 

which determines the actual value of <r; moreover we have 

sin* sin j, sin I /,^^ \ 

— • — = — '= ; . (109.) 

and it is easy to find <rco$s=p-\-p l cosI 1 (110.) 

<r cos i^Pj+p cos I J 

so that <r, sin/, sin/ y , cos/, cos/ ; are all simply expressible in terms of p, p f and I. 
The variation of a is easily found by means of the equation (105.), which gives 

^=(p / E i -^E)(Q / -Q).- (11 J.) 

The equation (103.) now gives (see (109.)) 

7 

<r(a/iCOSj>-- « sini')=^j(Q i — Q) ; . (U2.) 

and from (106.) we obtain 

<r 2 sinl(# o cosy+^isinj') = {(j0+^ . (113.) 

The two last equations determine the motion of the principal plane in space, irre- 
spectively of any arbitrary sliding motion which we may attribute to it in its own 
plane. For they give the angular velocities with which it is at any instant moving 
about two lines at right angles to one another in its own plane (see the end of art. 870- 
They may be put in another form as follows : — the actual value of Cl—CL is 

r r, 

j \ rpZ /yt/S 

where cos^ has the value given above (end of art. 86.); and if the operations indi- 
cated be actually performed, observing that Er=0, E y r=0, &c, and that 

Ecosx=— ^A &c, 
the results will be found to be 

(T T \ 

^ — 2 ) sin (0— v) sin (0, — p) 

(T T \ 

^i— -~2 ) X (p cos (0— p) sin (6,—v) +^ y sin (tf — p) cos {6 t — p)). 

Here 0—v 9 d t —p represent, it will be remembered, the angular distances of the planets 
from the line of nodes. 

We will assume for the present the condition &> 2 =.Q, so that the plane of xy may 
have no sliding motion, but roll upon the conical surface to which it is always a tan- 

* Referring to the arrangement supposed in the diagram, it will be seen that *■ becomes negative in the case 
now considered. 



mmA-z— isicosjc, 
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gent plane. If, then, we call at the actual angular velocity of the principal plane 

about its instantaneous axis of rotation, so that ^=a/^o+^i? and put J for the angle 
between the line of nodes and the instantaneous axis (which is the line in which the 
principal plane is intersected by its consecutive), we shall have 

# cos !>+*>! sin j>==a>cosj and ^cosv — co sm i>=a>sin/; 

and the above equations give 

ffCOtj=zpQOt(6—'V)+p l COt(0 l — v) 9 

a result which may also be put in the form (see (109.)) 

sinIcotj=sin/ / cot(^— *) — sin * cot (0,— v). 

It is very easy to show by spherical trigonometry that this equation signifies that 
the instantaneous axis coincides with the line in which the principal plane is inter- 
sected by the plane of the two radii vectores # . In other words, we have this theorem ; 

The principal plane always turns about the line in which it is intersected hy the plane 
of the two radii vectores. 

It follows of course that the principal plane always touches the conical surface 
described in fixed space (relatively to the sun) by the said line. I think it probable 
that most persons would expect, at first sight, that the principal plane would always 
touch the conical surface described by the line of nodes, which, as has been just 
shown, is not the case. It is perhaps worth while to verify this result by inde- 
pendent reasoning. 

89. Let Roman letters refer, as in former articles, to axes of coordinates having any 
arbitrary fixed directions. Then, putting A=m(yz'— zy'), &c. (as in art. 80.), and 
using §, i\, £ for current coordinates, the equation to the principal plane is 

(A+A^+(B+B>+(C+CX=0; (114.) 

and the line in which this plane is cut by its consecutive is determined by combining 
the equation (114.) with that obtained from it by differentiating the coefficients of 
§,.??, £ with respect to t\ namely, 

(a , +a;)?+(b / +b / >+(c , +cx=o (ii5.) 

Now from the fundamental equations 

.. x ail 

/ „ „ da da , wV _ * 

we obtain A'=ro(yz' — zjr^y-^- — z^==mm i (yz i --zy i )(8 3 — r~ 3 ) 

* Let $ be the angle between the principal plane and the plane of the radii vectores ; the former plane 
divides the angle I into two parts, of which one is i k and the other is — t; and we get two spherical triangles 
which have a common side/, with adjacent angles ^ and \{/ in one, — i and ir— -y in the other; hence 

cot (8 — y)sin/= cot rf/ sin ; + cos/ cos i 
cot(0,— -y)sin 4 /=cot vj/ sin ^ + cos/ cos i t ; 

and if cot \J> be eliminated between these, the result is the equation in the text. * 
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(where 5 is the distance between m and m,); and in like manner 

A'==— mmfyZt— zy,)(S- 3 — r" 3 ) ; 
we have therefore A'+A—mm^yz,— zy^r/ -3 — r~ 3 ), 

with similar expressions for B'+B|, C'+CJ; so that, when the common factor 
w/r^-r 3 ) is 'omitted, the equation (115.) becomes 

(y z /-~ z y^+( zx /- xz >+(xy / -yxX=0 ? 

which is evidently the equation to the plane containing the two radii vectores. Thus 
the theorem in question is verified. 

90. To return from this. digression : the motion of the line of nodes in the principal 
plane will be given by putting a 2 =0 in either of the equations (102.), and introducing 
the value of ^siny— a^cosv from (103.). In this way we find, after slight reductions, 

^siniys^cosi^+^cosi-^, 

in which we may substitute for cos/, cos/ ; , the values given by equations (110.); this 
gives 

cpp l sin I . *' = (p? +pp t cos I) ^j + (p 2 +PP,gosI)-^i 

or, if we introduce the actual values of ^p -^p we find 

pp>J= — mmprpm {6—v) sin (d l —p)x{p cos / . (cS~ 3 — r" 3 ) H-j^cos i$~ 3 — r~ 2 ) } . 

It is not my purpose however to enter further into details; and I shall conclude 
this subject by briefly examining the consequences of a slightly different assumption 
as to the motion of the axes of coordinates. I shall suppose, namely, that the plane 
of xy still always coincides with the principal plane, but has a sliding motion such 
that the axis of x always coincides with the line of nodes, 

91. The assumption made at the end of the last article implies the condition v=0 ; 
and oq 2 will no longer be 0, but must be determined by equations (102.) ; either of 
these gives (putting j/=/==0, and reducing by means of (103.), (109.), &c.) 

pp I smI.af 2 ==p i cosi i 'j^+pcosi-^, V 11 "-) 

which coincides, as of course it ought, with the expression given for J on the former 
hypothesis (art. 90.). The difference is that Q, Q, now no longer contain v. 

The values of a 09 ^ are obtained at once by putting ^=0 in the equations (112.), 
(113.) ; and all the conclusions which were derived independently of any supposition 
as to the value of <y 2 , subsist as before, when modified by putting v=0. 

We may add one more equation, which is required in forming some of the expres- 
sions for the variations of the elements ; namely, 

^ * — »|Sin I /n>7\ 

tan-= -' t* • * v 11 /*/ 

2 a- +p -\-pt cos I 
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This is easily obtained from (109.) and (110.) ; and in like manner 



tan 4 



^ p sin I 



2 cr -fp y -±-p cos I 

I shall now recapitulate the results of the last supposition, so as to exhibit in one 
view the transformed differential equations of the problem of three bodies. It will 
be as well to repeat also the explanation of the symbols. 

92. Signification of the Symbols. 
M, m, m t are the masses of the sun, and of the two planets, and M+mm^, 

a and e are the semiaxis and excentricity of the instantaneous ellipse described by 
m about the sun. 

a t and e { have the same meanings with reference to m r 

I is the inclination of the plane of the former ellipse to that of the latter. 

0, 6 l are the longitudes of the two planets, measured in the planes of their orbits 
from the common line of nodes. 

r, r t their radii vectores. 

t*r, w t the longitudes of the perihelia, measured likewise in the planes of the orbits 
from the line of nodes. 

5, Sj two elements such that I ndt-\-s, 4 n i dt-\-e l are the mean longitudes, where n 9 n t 

Jo Jo 

are defined as usual by the equations n 2 =^a"" 3 , wf=^a" s . 

Then 1 ndt+e— ns is the mean anomaly of m, and r, 6 are functions of the mean 

anomaly and mean longitude given by the laws of elliptic motion. The same is true 
for m t , mutatis mutandis. 

% is the angle between the radii vectores, so that 

cos %= cos 6 cos y -f sin 6 sin 6 t cos I. 

Let § be the distance between the planets, so that 

S 2 =r 2 +r 2 -- 2rr) cos %. 

Q, Q ; are the disturbing functions, defined by the equations 



,1 rcosy 
Q = mm,( -s:— ■ 



,_ /I r/CosvX 

, Q l =mm i (^^ J -^J, 



i® r 



& r 2 

and, when expressed in terms of the elements and t, are functions of 

r*t r*t ft nt 

a, a p e, e p I, 4 ndt-\-e, 4 ndt-\-e>--zF, 4 n i dt-\-e p 4 n^t-\-i i 

Jo Jo Jo 'Jo 

When Q, Q, are considered on the one hand as expressed in this way, and on the 
other, in their original form as functions of r, r„ 0, P and I, we have, as applied to 
either of them, 

d d d d d d 
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/?, p p c are defined by the equations 

/, / ; are the angles between the principal plane and the planes of the two orbits, and 
are given functions of the elements a, a p e, e p I, by virtue of the equations 

— per sin i=p i <r sin ^'=-pp l sin I, 
from which we have also 

(7 cos i~p-\-p i cos I, a cos i t =Pj+p cos I 
(for the values of tan ~? tanX see art. 91. ). 

a> is the angular velocity of the principal plane about the line of nodes ; 

a>! the angular velocity of the principal plane about a line in itself perpendicular to 
the line of nodes ; 

oo 2 the angular velocity of the line of nodes estimated in the direction of the prin- 
cipal plane. 



Differential Equations of the Problem. 

93. The nine intrinsic elements, as we may perhaps appropriately call them, namely, 

a, a p e, e p e, e p w, t*t p I, 

are determined as functions of t by the following system of nine simultaneous differ- 
ential equations of the first order : 

mpa* =2na 2 — •> m l pp=2n j a l ( -j- i 



m^e = 



mpe 



_nj?y^m +{1 _ vl 



[div i 



<) 



MO, 

de { 



2na^+ mVl - e \i _Vi Z^f 1 



da 



de 



sin I 



cos I dQ, 1 d£l { 
p dl^p t dl 



i ~ 2 dffl, , n,a, Vl—e 2 , *— 

m th s = -2ntf- s -l+JJ^ ,(] -^/j 

m, ft, f cos I d£l, 1 d£l 



e ' ] de, 



in 1 1 j» y dl *p dl 



sin 



m^ 



/ na Vl — e 2, d£l ma, 



de sin I 



cos I d£l , 1 rfll, 

L p dl l p t lT, 



_ / w.#, i/l— • e 2 iO, m.u,, [cos I d£l, , 1 dOl 
^ ' J de t sinll Pl dl^p dlj 
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PPi sm I.I'=^cos I(j + fa) + &;+3^; 



+p 



cos iffi + ^ + ^ + ^ 

I ^£ y * e/'STy / * de * dm 



94. The only parts' of the preceding expressions of which the deduction is not per- 
fectly obvious, are the terms involving I in the values of s r , s', sr', «r. They are 
obtained, as has been sufficiently explained, from the expressions in art. 55, to which 

are to be added the values of -rp &c. ((95.), art. 79-) ; °& putting v=0, a> disappears 

from the latter; and the values of o/ 1? a/ 2 , cos/, cos/,, tan^U tan^ are to be introduced 

(equations (110.), (112.), (116.), (117.))- After some rather troublesome reductions, 
the expressions above given will be found. 

In these equations it will be recollected that the mean longitude of m is represented 

C l d£l 

by I ndt-{~e, and the differentiation with respect to a in ^ is only performed so far 

as a appears explicitly. If we wished that the mean longitude should be expressed 
by nt+e, the only change in the equations would be that the differentiation with 
respect to a must be total ; i. e. must extend to a as contained in n, A similar 
remark applies of course to s r 

In actual use it would be more convenient to introduce R, R, instead of Q, Q,; the 
latter functions give a rather more symmetrical form to the equations, and are more 
convenient in general investigations. (The relation between them here is merely 
Q=mR, Q^w^R,; in another part of the paper the symbol Q was used for -mR 
(art. 55.)). 

95. If the equations of art. 93. were completely integrated, the intrinsic motion 
of the system would be completely determined; that is, we should know at any 
instant the dimensions of the two orbits, the mutual inclination of their planes, the 
position of their axes with respect to the line of nodes, the place of each planet in its 
orbit, and (by (110.)) the inclination of each orbit to the principal plane. 

The position of the system relatively to fixed space would then have to be sepa- 
rately determined as follows : — 

The three quantities cj , a> 19 a 2 (see end of art. 92.), of which the values are ((112.)* 
(113.), (116.)) given by the equations 

a sin 1.^=^+008 1)^+^+ cos i)^, 

MDCCCLV. 3 A 
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would be given functions of t. Then if we call 

J the inclination of the principal plane to an arbitrary fixed plane ; 

Q the longitude of the line of intersection of these two planes, reckoned in the 
fixed plane from a fixed line ; 

N the angle between this line of intersection and the line of nodes; 
we should have (as in art. 83. with a different notation) 

J'=fiJ COSN-- ^ sin N 1 

Q'sin J=fi/ sinN+a>iCosN I (118.) 

N'=fi; 2 — cot J(u sin N+^i cos N) J 

and the integration of this system would give J, Q, N as functions of £, and so deter- 
mine at any instant the position of the principal plane and of the line of nodes, rela- 
tively to fixed space. 

With respect to the motion of the principal plane, the following may be added. 
It has already been shown (art. 88, 89.) that the line about which it is at any instant 
turning, coincides with that in which it is intersected by the plane of the radii vectores ; 
and the values of a> , a x (see art. 88, putting »=0 in the expressions there given) may 
be put in the form 

a 2 w =mm l sin I . I ^— - 2 ) (p cos 6 sin 6 i -\-p i sin 6 cos 6 t ) 

trw l =zmm l &ml(-{ — - 2 ) sin Hmd r 

If the latter of these be multiplied by <r, and then both sides of each squared, and the 
results added (after putting for <r 2 on the right its value p 2 +p*+2pp l cosI), we find, 
observing that cos 6 cos t -\- sin 6 sin 6 i cos 1= cos '%, 

a 2 s/afl+^=zmm t sin I M^- ~ J (p 2 sin 2 0,+pf sin* 6+2pp t sin 6 sin 6 l cos x)*, 

an expression which may be further transformed as follows. Let X, \ be the latitudes 
ofm^^with reference to the principal plane) ; then sin X= sin 6 sin/, sinX^sin^sin/^; 

hence, since p= . ' » p,=- — r— ^— ? we obtain 

$s/o%+a\=mml-{ — ^J(sin 2 X+ sin 2 ^— 2sinXsinX / cos%)\ 

This gives the absolute angular velocity with which the principal plane is at any 
instant changing its direction in space; it is evident that (if the supposition r=r be 
excluded) it can never vanish except when both planets are in the line of nodes. The 
direction of the rotation is determined by the signs of <y and <y x . 

96. The system of differential equations given in art. 93. affords an example of the 
so-called "elimination of the nodes" in the problem of three bodies. Jacobi, by a 
very remarkable and ingenious transformation- has effected the elimination in a quite 
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different manner*. The equations of art. 93. are merely transformations of the 
original differential equations of the problem, without any integrations ; they are 
however in a form which might perhaps be used advantageously in certain cases for 
the purposes of physical astronomy. Those of Jacobi are obtained by employing all 
the four usual integrals of the problem, and are shown to include an additional inte- 
gration. They have however the disadvantage of substituting the coordinates of two 
fictitious bodies for those of the actual planets, and would probably be inconvenient 
for ordinary practical use ; though in a theoretical point of view they seem to deserve 
more attention than they have hitherto received. It would be wrong to take leave 
of this celebrated problem without referring to another transformation by M. Ber- 
TRAND-f-, which, as has been remarked by a recent writer in the same journal, effects 
six integrations, and therefore represents the furthest advance which has yet been 
made towards a rigorous solution. 

Appendix A. 

When the method described in Theorem VII. (art. 49.) is applied to the solution 
of a system of equations of the form (L), of which n integrals, a Y ... a n9 satisfying 
the conditions [a i9 %]=0, are given, the first step is to express the n+l partial dif- 
ferential coefficients ^- ? &c, and -^-; namely, the values of y u ...y n and —■ Z in 

terms of x u ... x n9 a l9 ... a n and t. The direct process is then to find X by integrating 
the expression y x dx x -\-y^dx % ^ ..-\-y n dx n — Zdt 9 and afterwards to form the remaining 

//IT 

integral equations ~^-=^b l9 &c. : when this process is adopted, the inferior limits in 

the integrations are perfectly arbitrary ; in other words, we may add to X an arbitrary 
function of a l9 a 29 ... a n9 without altering any of the general properties of the final 
system of integrals. 

But it is generally much more convenient to perform the differentiations with 
respect to a l9 ... a n% first 9 and integrate afterwards; thus we obtain the remaining 
equations in the form 



b-((&dx+&dx+ +&dx-—Jt\ 



When this plan is followed, the limits are still arbitrary if it be only required that 
the equations thus obtained shall be true; but if it be required that they shall form, 
with the given integrals, a normal solution, it is necessary to take the limits in such a 
manner that the functions equated to b l9 b 29 ... b n shall be the partial differential coef- 
ficients with respect to a l9 a 29 ... a n9 of one and the same function ; which will not 
generally be the case unless care be taken that it should be so. 

Iri practice, the expression for rfX usually consists of several terms, of which each 

* Comptes Rendus, 18.42, part 2. p. 236, &c. t Liouville's Journal, 1852. 

3 a 2 
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contains one of the variables only. Suppose one of these terms is 

VL/ltX, £*13 ^29 ••• UmJCI/JL^ 

so that, so far as this term is concerned, we have 



X=l <p(x 9 a 19 a 29 ... a n )dx. 

Ja 



where A is an arbitrary function of a l9 ... a n . Consequently 

and we see that we should not *w general obtain the differential coefficients with 
respect to a 19 ... a n of one and the same function X, by merely integrating --£? -^? &c, 

da y dcL<£ 

with respect to x 9 from the same inferior limit A, chosen at hazard. 

But it is evident that we shall attain this end if we adopt the following simple 
rule : — 

Integrate -=£■? &c. with respect to x 9 taking the same inferior limit in each ca$e, 

lA/Cl-i 

namely, either 

(1) a value A of x which satisfies the equation <p(x 9 a 19 ... a n )=0, or 

(2) any determinate constant (i.e. not a function of a 19 ... a n ). 

For example, in the problem of central forces (Part I., art. 28, &c), we had (see art. 19.) 

dX——hdt+cd0+(2m(h+<p(r)) — /fr-*y^ 

(where r, 0, X are the three variables). 

The very troublesome process of differentiating X with respect to h 9 k and c after 
first finding X by integrating the above expression, is avoided by the method adopted 
in art. 29 ; namely, by differentiating first, and integrating afterwards. In the inte- 
grations with respect to r, the inferior limit is one of the roots of the equation 

2m(h+<p(r)) — l?r-*=0, 

namely (in the case of elliptic motion), the perihelion distance; and in those with 
respect to X 9 the inferior limit is ; so that the rule above given is observed. 

At the time of writing the article referred to, neither the rule itself, nor the neces- 
sity of attending to the limits, had occurred to me ; it was therefore, strictly speak- 
ing, accidental that the final integrals were obtained in a normal form. 

In treating the problem of rotation (Section III.), I perceived the necessity of cau- 
tion as to the limits, if the former order of proceeding were adopted ; but preferred 
avoiding the risk of error altogether, by performing the integrations first, so as to 
obtain the actual expression for V. The final equations (R.) 9 art. 44, might however 
be obtained in a more simple way by differentiatingj^raf ; thus we should have (see 

77 7 

equations (45.), (46.)), observing that--=— , &c. (art. 44.), and putting 

Ullh Ziflr 
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(1— cos 2 /— cos 2 j+2cosicosjcos#— cos 2 #)^=Q, 
dV k 



dh~~2h 
dV 



k 



dco&j' 



4* cos i+(p oosj+\-^^dd\=t-\-r 
f . r cos^ cos fl— cosi 

JVT] — 

f+ Jsin0 



dcosi 1 J Qsin0 

dV {C-A)k 3 co$j 



(i.) 



2ACA 



^cosi+<pcosj 



^ 



• 7 f i T cos ? * cos ^ — cos i 7/i1 ■ /> 

^j^» we 

may (1) eliminate this term between the first and last of these equations, and (2) 
eliminate 4 between the first and second. We thus find the two following equations, 



, Tcosicos^— cos? 7A & C-— A, ./,, x 



(ii.) 



cosjl <p + 1 



cos i cos 0— cos^* 
QsinS 



dA 



acosz 



. Tsin dA 2h,, f x ~ _. „ 



which last, combined with the preceding, gives 



rsin 6dS k, x a cos i + j3 cosj" 
1 ~Q~ == A^+ r ) 






A 



(iii.) 



and we may take (i.), (ii.) and (iii.) as expressing the solution of the problem. 

Now we have cos I=— Z "T cos l^- } from which it is easy to find (observing the 

sinj sin o ■ J v ° 

conditions which determine the sign of Q) 

7T cos i cos 0-— cos; 1A t . ., , 

— dl= tt-^—t—* do, and similarly, 

H sm v 9 J> 

j T cos / cos — cos i j 
aJ = ^ .♦ . A dv ; 



Qsin0 



we have moreover 



d0= 



sin 6d6 



a 



All the integrations may therefore now be performed immediately; and we may 
take for the inferior limit of 6 any value which satisfies the equation Q=0, or 

(cos 0— cos i cosj) 2 — sin 2 i sin 2 j=0 ; 

this is satisfied by 0=i-\-j 9 which evidently corresponds to 1=0, J=0, 0=0, and it 
is manifest that equations (i.), (ii.), (iii.) will thus become identical with equations 
(R) of art. 44. 

I do not regret however having introduced the rather prolix investigation of arts. 
39 and 43, because it is interesting to know the actual value of V (equation (48.)), 
which the method just given leaves undetermined. 
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Appendix B. 

On the subject of the transformation of elements, the following additional remarks 
will hardly be superfluous. Suppose Q is originally a function of the elements a, 5, c, 
&c, with t; and let m, |3, y, &c. be other quantities connected with a y b, c, &c., by 
equations such as 

da=Adt*+Bdp+Cdy+...+Kdt 9 • . (a.) 

r 

where A, B, C, ... K are given functions of a, |3, y, .... £; or by equations such as 

rfa=A i £fo+B / rf6+...+K / rf^ ........ (a.) 

where A /? &c. are given functions of a, b, ..,., t. In either case, if each of the equa- 
tions be integrable per se, we may consider a, b, c, .... as functions of a, (3, y, ...., £; 
and such equations as 

are both significant and true. 

But if the expressions on the right of the equations (a.) be not differentials per se, 

the equations (Q.) are either unmeaning or untrue. For the symbol ~~r implies one 

of two things; either, that Q is expressed in terms of a, j3, ... £ without arbitrary 
constants (i. e. that the transformation of Q can be actually effected without inte- 
grating the differential equations of the problem), which is manifestly impossible, 
unless (a.), &c. be integrable per se; or else, that the differential equations are to be 
conceived to have been completely integrated, so that a, i, &a, and consequently 
A p B /? &c, are known as functions of t and arbitrary constants, whereby the right- 
hand side of (a.) becomes an explicit function oft (and arbitrary constants), so that 
&, j3, &c. may by integration be expressed in the same way, and, by means of (a.), 
a, b, &c. may be similarly expressed, and finally, by algebraical elimination, a, b, &c. 

become functions of a, /3, &c, t, and arbitrary constants. On this supposition, ^ 
has a meaning, but the equation (O.) is untrue; for we must have 

da"~" da daT db ^+~** ; 
da 

and it is manifestly not true that ^=A, &c. in this case, because the equation 

da=Adm+Bd(Z+...+Kdt, 

not being integrable per se, only subsists for those variations of ^, |3, &c. which 
actually take place during the instant dt ; whereas the equation 

A 

subsists for arbitrary variations of all the variables. This view of the subject entirely 
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coincides, in substance, with that taken by Jacobi ; but the above mode of stating it 
may tend to make it clearer, and to call attention to a matter which, so far as I know, 
is not so much as mentioned in any of the elementary works usually in the hands of 
students of physical astronomy. 

[Addition to Appendix B.] 

Received March 1, 1855. 

The remark made above, that the symbol — is unmeaning in the case considered, 

is not of course intended to imply that a meaning may not be given to it; but then 
such meaning is different from the ordinary signification of the symbol, which is a 
partial derived Junction. 

The whole matter may be strikingly illustrated by a simple example. 

Consider the movement of a rigid body about a fixed point. Adopting the nota- 
tion of art. 40 (Part I.), we have 

pdt= — cos <pdA— sin 6 sin <pd^ 

qdt= sin <pd6— sin 6 cos <pd^> 

rdt=d(p-{- cos 6d^. 

Let a, |3, y be three new variables, defined by the equations 



a 



nt Pt £t 

-\ pdt, j8 = i qdt, y=j rdt ; 

Jo Jo Jo 



so that du~pdt, &c. Then the above equations give 

dd= — cos <pda-\- sin <pdfi, rf<p=c?y+ cot 0(sin (pdu-{- cos <pdfi) 

rf-4/=: — cosec 0(sin <pdcc-{- cos tpdfi). 

Here #, j3, y are the sums of the elementary angles described about the axes in the 
course of the motion ; and no one would maintain that 6, <p, $ are functions of &, (3, y, 
for the values of the latter variables at any time do not determine the values of the 
former. If therefore we choose to write such equations as 

dQ dQ . D 

we must admit that -r> -^ &c. are not partial derived functions in the ordinary sense. 

da dj3 r J 

At most, -j- is the derived function of that function of a which 6 would become if /3 

and y were maintained invariable, i. e. if the motion were restricted to a rotation 

J 2 dH 
about the A-axis. Again, if we admit such symbols as -— -> -7-753 we must inter- 

pret them as follows : 

dH d dQ dcos<p_^ . d<p 

~dfiot = d^ d*~ = ~~ ~jf~— sm $dip 
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but jJ|= cot 6 cos (p, and therefore 



dH A . 

: cot 6 sin <p cos 9 ; 



and in like manner we find the same value for j-j# so that in this particular case 
the condition j^='jq^ * s verified. 

But if we take j-^ and ^rj- in the same way, we find the former = sin <p anrf £&e? 

latter =0, so that in this case the condition ' w not verified. The geometrical meaning 
of this is obvious ; analytically it is merely an instance of a general fact, pointed out 
by Jacobi ; namely, that the effect of two successive pseudo-differentiations with 
respect *to two independent variables, is not generally independent of the order of 
operations. 

If V be the potential of another body, given in position, upon the body considered, 
then V is a function of 6, <p, ^ and 

7tt dY 1 A . dV 7 . dV , ,. 

and if we substitute for dd, d<p y dip their values in terms of doc, ^/3, dy, we obtain an 
expression which we may call Lrfa+Mrf/3+Nrfy , L, M, N being functions of 0, <p, %£, 

~i6 9 ~r > -jt> °f which, as is well known, the mechanical meanings are the moments of 

ad d<p dy 

the attraction of the second body about the three axes. Here again no one would 

maintain that V is a function of a, |8 5 y ; and if, as is often done, we say -r- = L, &c, 

the above remarks apply in all respects to these equations. 

I should have thought it superfluous to dwell so much on these points if it had not 
appeared that writers on physical astronomy have in some instances either overlooked 
the distinction between real and psrarfo-drfferentiation, or at least have failed to point 
it out to their readers. The only discussion of the subject which I have met with is 
that given by Jacobi, in the correspondence referred to. 

It may be added, that in general investigations, where symbols such as -r- 9 &c. 

may be used without defining the nature of V, or the precise meaning of a, /3, &c, 

serious errors might be committed if it were assumed that the condition -r~iz=inrr 

dadp dpda 

always subsisted. 

Appendix C. 

The theorems relating to the transformation of coordinates, given in Section VI., 
may be made more general, and in many cases more useful, as follows : — 
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If x l9 x^ ... x n be the coordinates employed in the first statement of any dynamical 
problem, the differential equations are comprehended in the formula 

[If there be any forces, such as those arising* from a resisting medium, which do 
not satisfy the natural conditions of integrability, then on the right-hand side of the 
formula (D.), instead of we shall have an expression such as S^(X^) ; but such 
terms are easily introduced and allowed for separately, and do not affect the follow- 
ing investigation. I shall therefore here assume that they do not exist.] 

In the above formula, W is a given function of x l9 ... x n9 x l3 ... x n9 which may also 
explicitly contain t. 

In Section VI. the only case contemplated was that in which x 19 ... x n are inde- 
pendent coordinates ; in which case the formula (D.) is equivalent to n separate equa- 
tions, since $x 19 &c. are wholly arbitrary and independent. 

In practice., however, it is often more convenient to use, at first, a set of coordinates 
more in number than the independent variables of the problem, and therefore subject 
to equations of condition. 

Let us assume then that the n coordinates x l9 ... x n are subject to r equations of 

condition, 

Lj=0, L 2 =0, ... L r =0, 

where L x , &c. may explicitly contain t 9 besides the n variables x 19 &c. 

dW 
If we introduce the n conjugate variables y l9 ... y n defined by the equations y^-j-r? 

and take Z a function of x u &c, y 19 &c. (with or without t), defined by the equation 

Z= — [W] + [x l ~]y l + . . . + [oQjjn 

(the brackets indicating that«ri, &c. are expressed in terms oiy l9 &c), then it follows 
exactly as in art. 18 (Part I.), that the formula (D.) will be changed into the system 

dZ 






dy { 



^(^+^|)^= 



(E.) 



[In the most usual problems W is of the form T+U, where T is homogeneous and 
of the second degree in x l9 &c, and U does not contain x l9 &c. at all. In this case 
Z is only T—U expressed in terms of -y l9 &c., instead of x l9 &c. But T is not neces- 
sarily homogeneous ; in fact it is not so in problems relating to motion relative to 
the earth, as affected by the earth's rotation.] 

Let us now suppose that the system (E.) is to be transformed by the introduction 

of the m independent coordinates | 1? f 2 , ... % m9 and of the new conjugate variables 

^15^25 — n m \ where m=n—r. And let it be required to investigate a theorem by 

means of which the transformation may be effected without recurring to the original 

formula (D.). 

MDCCCLV. 3 B 
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The definitions of the new coordinates g„ &c. will furnish m equations (which may 
explicitly contain t) by means of which | 13 ... l m may be expressed as functions of 
x u ... x n (with or without t) ; and conversely, by means of these m equations, together 
with the r equations of condition 1^=0, &c, the n variables x l9 x 2 , ... x n may be 
expressed as functions of £„ ... | m5 with or without t. When x 19 ... x n are so expressed^ 
let them be represented by (a?,), ... (x n ). We shall have then 

l _d{x l ) d{xj) , .d(%i) _, f 

Xi ~~~dt "+" d% x * 1 "^''' ' C* m ' •..•••• \X.) 

so that x[, &c. are expressible (and in only one way) as functions of fj„ &c, fi, &c. 

If then the formula (D.) be transformed by expressing^, &c, x u &c. in this manner, 
it becomes, as is well known, 



-m- d wh=°' 



where (W) represents the result of transforming W as above ; and since $% 19 &c. are 
now independent, this formula breaks up into the m separate equations 

w(W)\' d(W) 



d(W) 
Moreover, if we now define % by the equation ~ } ^r l =7i i , and put 



. . . (F.) 



where (Si), &c. are expressed in terms of n l9 &c, we know already (art. 18.) that the 
system (F.) becomes 

^ i= ~(fa' ni ~~'Wi ' ' 

Now let P be a function of the m new variables | 15 ... % m , and of the w old variables 
#i> ••• y w (with or without t) 9 defined by the equation 

P = [x 1 )y x + {x 2 )y 2 + . . . + (x n )y n . 

d(W 
Since ^=n-^r 9 and since (W) contains ii, &c, only through x' i9 &c, we have, ob- 



. ^W 
serving that -^r—y i9 



«() — — flf f I 7/ . — 4— r L >J/ ** 



but from the equation (a?'.) we have ^==.-~^ 

dix 1 dix ) dix ) 

consequently * = #» d&" +y»~ifif +••'+&« dt ! 

dP 

an expression evidently equivalent to ^- Thus ^ mai/ &e defined by the equation 

— — 

jy m ~—-*li • • • * » » » » « » * » » » e \fj») 
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without recurring to the formula (D.). And since each equation of condition gives, 

dZ 
if differentiated totally with respect to t, with the substitution of j- for x i9 an equa- 
tion such as 

OdLt dlu dZ ttL dZ , T 

— dt^dx l dyC r ~^ r dx n dy n v^v 

the m equations (*;.) with the r equations L (in which last x l9 &c. must be expressed, 
after the differentiation, in terms of |, &c), give n equations by means of which 
y 19 ... y n may be expressed in terms of | 15 &c, jj 19 &c, and can be so expressed only in 
one ivay. 

Lastly, the value of Y (see equations (G.)), may be obtained as follows: — • 

Since Y=-(W) + (? I )*.+ ... + (£.H. 

and Z = — [ W] + \jc\~] #!+...+ [#»]y« 

and (W) is only [W] differently expressed, we have, without reference to modes of 
expression, 

t — Zj "=• z>i [ %$i) —■ ^i \ $iyi) • 
On the other hand, since P can contain t explicitly only through (x x ) P &c, we have 

dt—^yi dt y 

and also H^)=^r§) + 2 *(ff 

dP 

hence, observing that ^=^ ? we obtain 

dV 

2«(g%) — 2,(%) = — -# * 

and therefore, finally, "¥=Z—--?t> 

so that Z-~-fo must become identical with "¥*, when expressed entirely in terms of 

the new variables. 

These results may be stated in the form of the following theorem. The system (E.) 
is transformed into the system (G.) by the following substitutions : — 

(1) a?„ ... x n are expressed in terms of f 1? ... l m by means of the m equations which 
define the latter variables, together with the n—m equations of condition 

Lj=0, ... Lv=0 (r=n— mi). 

dP 

(2) ij 19 ... jj m are defined by the m equations jjt==^; where the modulus of trans- 
formation P is given by the equation 

p =(^i)yi+(^)ya+—+Wyi., 

(a?,), &c. being expressed in terms of £„ &t\, so that P is explicitly a function of l l9 ... 
y\ "'Vm with or without t. 

2 B 2 
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dP 

(3) If is defined by the equation ' V F=Z— -p in which (after the explicit differentia- 
tion of P with respect to t) 9 #„ &c, y 19 &c. are to be expressed in terms of the new 

dP 
variables. y l9 &c. are thus expressible by the help of the m equations -717=^ and the 

dh /dL dZ\ ■ 
n-m equations ^+2^^j==0. 

dP 

If (^j), &e., do not contain t explicitly, then ~jt=0, and Y is obtained merely by 
expressing Z in terras of the new variables. 



It may be observed that the whole of the above reasoning would apply to the case 
in which the new variables | 1? ... % m are more in number than the independent vari- 

dP 
ables of the problem (or m>w— r), with this exception ; that the m equations ^f=^ ? 

together with the r equations obtained by differentiating the equations of condition 
totally with respect to t, would be more than sufficient to express y l9 ... y n in terms of 
the new variables ; consequently y l9 &c. might be so expressed in different ways, and 
therefore, although the value of Y obtained by the above rule would certainly be the 
same as that obtained by recurring to the original formula (D.), the form of If might 
be different, and therefore the resulting formula erroneous. 

There must doubtless exist some rule for choosing n — m combinations of the 
equations of condition in such away as to lead to the correct forms of y 19 ...y n as 
functions of the new variables ; but I have not at present attempted to investigate it, 
and perhaps it would be hardly worth while. The theorem in the case in which the 
new coordinates are independent, may, I believe, be practically useful. 



ERRATA IN PART I. 

Art. 1. equation (4.), for doc read dx{. 

Art. 10. In paragraph preceding equation (26.) omit the words " not containing t explicitly. " 

Art. 18. equation (/3), for yi read y { . 

Art. 19. equation (29.), for h read 6j. 

Art. 24. second line after equation (L.),/or " such as h, k" read " such as/, g" 

Art. 30. The expressions equated to h, k, c, and the three terms in the left-hand column of the table of 

elements, should each be multiplied by m. 
Art. 42. near the end, for " according as 6 is between o and ir, or not" read " according as is between 

7r and 27r, or between o and W." . 



